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SYNOPSIS 

"ANALYTIC TECHNIQUE FOR HYDRODYi'JAMIC INSTABILITIES 
WITH REALISTIC BOUTvIDARY CONDITIONS" 

A Thesis Submitted 

In Partial Fulfilment of the Requirements 
For the Degree of 

DOCTOR OF PHILOSOPHY 
By 

KRISHNA KIMAR 
to the 

Department of Physics 
Indian Institute of Technology, Kanpur 

September, 1986 

Hydrodynamic instabilities v/ith realistic boundaries 
are studied analytically in this dissertation. This is 
necessary for any~quantitative comparison between theory and 
experiment, the latter having become increasingly accurate 
over the last two decades. The effects of realistic boun- 
daries can be large for unmodulated flows and subtle for 
modulated ones. Flows are treated in Rayleighr-B-5nard (RB) 
and Couett e-Taylor (CT) geometries. Convective instabili- 
ties are discussed in both single component and binary fluids 
in Rayleigh-Benard configuration. 

The traditional boundary conditions assume free bounding 
surfaces, where the normal caaponent of velocity vanishes 
and there are no tangential stresses on the fluid. These 
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constraints make the study of hydrodynamic equations simple. 
However, these boundary conditions cannot be realized in 
practice. Realistic boundaries are, in general, rigid and 
impermeable. At rigid surfaces, not only the normal component, 
but also the tangential components of velocity vanish because 
of "no-slip” condition. The impermeability of realistic 
boundary surfaces prevents any mass current across them. 

These constraints make the analysis of hydrodynamic systems 
very difficult, A procedure, due to Chandrasekhar, which 
makes analytic study with rather general boundary conditions 
feasible, is followed here, 

Chandrasekhar's technique consists of Fourier expanding 
the field (usually the temperature fluctuation), on which the 
only requirement is the vanishing at the boundary surfaces, 
and truncating the expansion at a finite number of teirms. 

With this ansatz on the temperature fluctt^ation the remaining 
fields are exactly determined from the differential equations 
satisfied by them, and the critical Rayleigh number is evalu- 
ated self-consis tently from the equation for temperature 
fluctuation. In case of CT flow, the azimuthal velocity 
fluctuation is expanded and the critical Taylor number is 
determined from the equation for this fluctuation. The 
method yields accurate results in the lowest order truncation. 
The accuracy can be further improved by considering higher 
order terms of the expansion. In Chapter I of this thesis, 
a review of the procedure is presented by studying single 
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component fluid in RB geometry 

In Chapter II, the most general binary liquid mixture 

3 4 

( He- He) is studied in RB configuration using this techni- 
que, The critical Rayleigh number for the onset of station- 
ary convection is determined in the limit k/s» 1, v\tiere 
' k’ is the separation parameter and 's' the ratio of mass 
and thermal dif fusivities. The result is within seven 
percent of the exact numerical value in the lowest order 
truncation. The truncation of the expansion of temperature 
fluctuation at the next higher order reduces the error to 
one percent, While the traditional (free-free) boundary 
conditions lead to a stationary convection with a non— zero 
wave number, the realistic boundary conditions cause a 
stationary convection vdth zero wave number over a significant 
range of parameter values. 

The procedure is applied to investigate the hydro- 
magnetic convection in RB geometry in Chapter III, The possi- 
bility of oscillatory convection is explored in the limit of 
large Chandrasekhar number (Q >>1) and small magnetic Prandtl 
number (P2<<1). It is shown that the Cha ndrasekhar-Gibson 
criterion is not a necessary condition for overstability, 
but a necessary condition for the existence of a polycritical 
point , 

The onset of the first instability in a modulated 
single fluid in RB configuration is discussed in Chapter IV. 

The modulation of temperature difference between the two 
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plates delays the onset of convection, i.e., the correction 
to the Rayleigh number is positive. At high frequencies, 
the correction goes to zero as the fifth power of the frequ- 
ency of modulation* 

In Chapter V, this method is used to study the onset 
of convection in modulated binary liquids. The effect of 
modulation can be either stabilizing or destabilizing 
depending upon the pamneter values. 

The reported analytic and experimental results for 
modulated CT flow do not agree over the entire range of the 
frequency of modulation. Hall’s analytic study shows 
destabilization over the whole spectrum of frequencies. 

In Chapter VI, the modulated CT flow is studied using 
Chandrasekhar's procedure. The flow is destabilized in 
the low frequency regime and the results are in good agree- 
ment with that of Hall, In the very high frequency regime, 
the flow is stabilized a result, qualitatively different 
from Hall's theoretical 'work. The correction to Taylor 
number goes to zero as the fifth power of the frequency of 
modulation. The procedure, used in the thesis, thus also 
presents a systematic method for treating the effect of 
modulation in presence of realistic boundaries. 

In the appendix, the onset of oscillatory convection 
and subsequent period doubling bifurcations in modulated 
binary liquids are discussed using stroboscopic maps for 
the system. 



OVERVIEW 


The onset of hydrodynamic instability in fluids has 

1 “4 

been a subject of extensive study , The recent spate of 

5-7 

sophisticated experiments has focussed considerable 

8 

attention on externally modulated fluids . Analytic 
studies of the modulated systems are limited largely to the 
case of idealized boundaries, which are physically unrealis- 
tic. The results of a recent experiment reveals that even 
the qualitative features predicted by theory (mth idealized 
boundaries) are not all true for a hydrodynamic system vdth 
realistic boundaries , Thus for any quantitative comparison 
betweea theory and experiment, the hydrodynamic systems must 
be analyzed with realistic boundaries. In this dissertation 
the onset of hydr ociynamic instability in unmodulated as well 
as modulated flows vdth realistic boundaries is analytically 
studied in Rayleigh-Benard and Couette-Taylor geometries. 

Thermal convection in a single-component fluid in 

Rayleigh— Bilnard configuration represents the simplest 

9 

example of convective instability in fluids. In Rayleigh- 
Benard geometry, a thin fluid layer confined betvwsen two 
infinitely extended horizontal plates is heated fr cm below. 
As the temperature of the fluid is non— uniform throughout 
its volume, heat energy is conducted frcm hotter regions to 
cooler regions and a temoerature gradient is maintained 
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across the fluid layer. On account of positive 

thermal expansion, hotter parts of the fluid expand and 
become lighter than cooler parts. The arrangement being 
'top-heavy’ is unstable^ and the temperature gradient is 
called adverse. The fluid tends to redistribute itself 
to make the configuration stable (bottom-heavy). This 
natural tendency is opposed by the viscosity of the fluid; 

As the adverse temperature gradient is raised by raising 
the temperature of lower plate, the arrangement becomes more 
unstable^and for a critical value of the temperature gradient 
the fluid starts circulating. This causes the onset of 
thermal convection— ^n example of hydrodynamic instability. 

In the case of a single-component fluid, the velocity and 
temperature fields do not vary with time at the onset of 
instability; that is, the convection is stationary. 

By controlling the temperature difference between 
the two plates, the adverse temperature gradient can be 
controlled. So, the temperature difference between the two 
plates is the 'control parameter' of the problem. The 
response of the system to a time-periodic disturbance in 
the control parameter is important. It can lead to drastic 
modifications in the flow pattern. By varying the temper- 

I 

ature of the lov\er plate periodically, a time-periodic distur- | 
bance (an external modulation) can be imposed on the control | 

parameter. The external modulation of the temperature { 

’ t 

difference between the plates is known to delay the onset of | 



convection. This means that the external modulation raises 
the critical value of temperature difference and hence 
stabilizes the flo'.A^- pattern [Chapter IV], 

The critical temperature difference depends on the 
types of bounding surfaces considered. The use of traditi- 
onal boundary conditions, which assume free bounding surfaces, 
simplifies the mathematical complexities of theoretical 
studies, Hovyever, for any quantitative comparison with 
experiments, an analytic study must consider the fluid layer 

confined between realistic boundaries. Most of the theore- 

10 11 

tical studies with rigid boundaries use either Galerkin * 

12 1 
or Vi/.K.B. technique, Chandrasekhar's techniques 

(discussed in Chapter I) makes analytic study with rather 

general boundary conditions feasible, n ' pjiiwiiiNitiy 

This procedure yields 

quite accurate results for various convective as well as 
centrifugal instabilities in fluids, 

1 3—20 

A binary fluid (a mixture of two nonreactive 

miscible fluids) or a fluid with a solute (e,g, saline 
water) is a highly versatile hydrodynamic system for the 
study of convective instability. The tvyo-component system 
is characterized by two different diff usivities — " usually 
those of heat and mass (of the ^soMte ) — and, hence, also 

known as a double-diffusive system, A simultaneous presence 
of two different diffusion coefficients can 

i) drive convection even if the system is hydro- 
statistically stable^^ *^^^ and 
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ii) induce either stationary or time-periodic 

23 

(oscillatory) convection at the onset of 
of the first instability. 

The onset of convection in hydrostatically stable 

configuration cannot be realized in a single-component 

21 

fluid. This effect 'ms first anticipated by Stem 

22 21 

following an observation by Stommel et al, , Stern 
considered a thin layer of saline water with satit concentration 
increasing upwards heated from above , the configuration 
is also known as 'fingering regime’. The adverse concen- 
tration gradient tries to destabilize the arrangement, 
while the favourable temperature gradient tries to stabilize 
it. For a temperature high enough at the top, the density 
of the solution decreases upwards making the system hydro- 
statically stable. Now, if a small parcel of the solution 
is displaced upwards, it gains heat but not appreciable 
amount of salt from its surroundings because of small solute 
diffusivity compared to its heat diffusivity. The parcel 
becomes hotter and, therefore, lighter. It keeps moving 
upwards and convective instability occurs in the system. 

At the onset, convection is always stationary in this 
configuration. 

An oscillatory convection at the onset of hydrodynamic 

23 

instability was first realized by Veronis • In Veronis 
configuration also known as ’diffusing regime' a thin 
layer of saline water with its concentration decreasing 
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upwards is heated from below. If a small parcel of the 
solution is displaced upwards, it loses heat but nqt salt. 

It becomes heavier as it cools, and the buoyancy force 
drives it back towards its initial position. The parcel 
starts oscillating, A lag in temperature between the parcel 
and its surroundings makes the amplitude of oscillation grow 
in time provided viscosity of the solution is not too large. 

This means that the parcel returns to its mean position of 
rest faster than it leaves. The oscillatory instability 
thus produced is called overstable convection. The insta- 
bility is generally oscillating for saline water in this 
configuration. For other solutions, the first instability 
can be stationary also. These effects can occur also in a 

binary fluid mixt’ure , It was first demonstrated by Schechter 
14 

et al, . In fact, a system of binary liquids is more 

advantageous than salt-water mixture for experimental studies. 

3 4 

In case of binary liquids, particularly He- He mixture, the 
ratio of the two diff usivities can be parametrically varied 
over a vd.de range by varying the temperature difference between 
the plates. It is very d3,fficult to do the same experimen- 
tally for a salt-water mixture. 

Most of the theoretical studies on binary liquids are 
confined to idealized (free) boundary conditions , The effects 
of realistic (rigid and impermeable) boundary conditions can 
be drastic. While the free boundaries cause stationary convection 
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with non— zero wave number, realistic boundaries lead to 
stationary convection of zero wave number over a significant 
range of parameter values [Chapter II j. If an external 
modulation is imposed on the temperature d 5 .fference between 
the bounding surfaces, the onset of stationary convection can 
either be delayed or preponed depending upon the parametric 
values. That is, the modulation can either stabilize or desta- 
bilize the flovi; pattern in binary fluids depending upon the 
parametric range. In case of oscillatory convection in binary 
fluids, the modulation a Ivways destabilizes the flow— pattern 
[Chapter v]* 

If a single-component fluid enclosed in Rayleigh-Benard 
configuration is electrically conducting, an externally imposed 
magnetic field delays the onset of thermal instability. At 
the onset the convection is generally stationary. The possibi- 
lity of a time-periodic (oscillatory) convection in this system 

1 

is a problem of long standing. Chandrasekhar studied this 
problem using idealized boundary conditions and showed that 

the magnetic Prandtl number ( P2) hsd to be larger than the 

thermal Prandtl number (p^) for any possibility of oscillatory 

convection at all. However^ under terrestrial conditions 

24 

P2 <<P-|. Later, Gibson analyzed this hydromagnetic system 

with realistic boundary conditions and obtained Chandrasekhar's 

25 

criterion once again. Recently Banerjee et,al, predicted 
the possibility of overstable convection — oscillatory 
convection with growing amplitude, for P2 < P-j 
but in the presence of large magnetic field. This 
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demands re-examination of the flow with proper boundary 
conditions. In Chapter III, it is shown that Chandrasekhar- 
Gibson criterion is not a necessary condition for the onset 
of oscillatory condition, but rathex a naceaaary ’crondi'tion for 
the existence of a polycritical point — ^ a point in para- 
meter space where the line of oscillatory convection meets 
the line of stationary convection, 

Couette— Taylor flow is the simplest system to study 

the hydrodynamic instability in rotating fluids. In this 

arrangement, a fluid is enclosed between two infinitely long 

coaxial vertical cylinders. The inner one of them is 

rotating with a constant angular velocity (Q), w^iile the 

outer one is at rest. Since the inner layers are rotating 

faster than the outer ones, the fluid tends to move radially < 

This is counteracted by a pressure gradient developed because 

of viscous drag. There is no vertical component of the 

velocity for lo’w rotational speed because of the symmetry of 

the system. Thus the basic flow, also known as ’Couette flow’^ 

is always tangential to the cylinders. As ’Q’ is increased, 

the centrifugal forces on inner layers grow, and for a critical 

^ ( = they become larger than the viscosity induced pressure 
^ % 

gradient. As the inner layers cannot move outward uniformly 
because of the crutter layer in their ways, they break into 
bands and circulation starts. These bands are called 'Taylor 
Vortices' and the transition from Couette flow is termed as 
’ Couett e-Taylor instability'. At the onset of instability. 
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the flow pattern is stationary. Further enhancement of 'Q' 
makes Taylor vortices wavy and finally leads to turbulence 
in the fluid. 


The reported theoretical and experimental results 

26-30 

for modulated Couette— Taylor flow do not agree over the 

entire range of the frequency of modulation. The full 
hydrodynamic equations for this system have been studied 

26 2n 

analytically by Hall and numerically by Riley and Laurence, 

IflJhile Hall finds destabilization of flow due to modulation 

29 

at all frequencies, Riley and Laurence notice that, for 

large amplitudes of modulation, the effect of modulation 

is to destabilize the flow for low frequencies and stabilize 

it for high frequencies. For small amplitudes, Riley and 

Laurence do not notice the stabilization at high frequencies 

perhaps because the effect at high frequencies is very small 

and barely within the range of numerical accuracy for moderate 

amplitude and beyond it for smaller amplitudes. The 

29 

experiments of Donnelly show a stabilization at all frequen- 

30 

cies, while the more accurate experiments of Thompson show 
a destabilization at low frequencies and a stabilization at 
high frequencies. In Chapter VI of the thesis we analyze 
the problem with realistic boundary conditions, using 
Chandrasekhar’s technique, and show that the effect of 
modulation (small amplitude modulation) is to destabilize 
the flow at low frequencies and to stabilize at high 
frequencies a result qualitativesiiy different from that 
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26 

of Hall . 

31 —34 

A low-order truncation facilitates the handling 

of non-linearity in hydrodynamic equations. In the appendix, 

35 

a stroboscopic map is constructed from Lorenz-like trun- 
cated equations for binary liquids and the possibility of 

36-37 

period-doubling chaos is discussed. 
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CHAPTER I 


A REVIEW OF CHANDRASEKHAR’S TECHr'IIQUE 

1,1 Introduction 

1—3 

The advent of sophisticated experimental technicfues 

has provoked considerable attention recently in the study of 

4-7 2 

hydrodynamic instabilities • Recent experimental results 

show that even the qualitative features predicted by theory 
are not all shared by the hydrodynamic system mth realistic 
boundaries. Thus to obtain any reliable results, the effects 
of realistic boundary conditions must be included in the 
theoretical studies of hydrodynamic instability. The consider- 
ation of "no-slip" condition because of the rigidity of 
realistic boundaries makes the analysis very difficult, 
particularly in presence of an external modulation. 

The effect of realistic boundaries on both modulated 

and unmodulated hydrodynamic flows cam be studied by using 
4 

Chandrasekhar's technique. In this chapter, we review the 
technique by studying the convective instability in single- 
component fluid in Rayleigh— Ben a rd geometry. The instability 
in the modulated flow has been studied in Chapter IV, 
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1,2 Hydrodynamic Problem :• 

An iacompressible fluid of viscosity ’p' is confined 

% 

between two rigid horizontal plates of infinite extension. 

The plates are a distance 'd’ apart. The lower plate is 
at a temperature T^ , while the upper one is at T 2 ( < T.^ ) , 

In conduction state there is no motion in the fluid. This 
means that the force due to gravity is balanced by the 
pressure gradient,- Steady state pressure field is determined 
from the equation 



P g 
s 


(1 .1 ) 


and the conduction of heat in the same stats is expressed as 

2 

3 T^ 

■—4“ = 0* (1 -2) 

3 z 

The boundary conditions on temperature field determines T , 
vjhich found to be 


Ts = T (-f) , (1 .3) 

where T^^ (= is the mean temperature, AT (= T^-T 2 )» 

the temperature difference between the plates. The 
co-ordinate system is such that the plates are normal to • 
z-axis at positions z = +1/2 and z = -1/2, 

As v>;e raise the temperature difference (A T) , for a 
critical value convection starts , At the onset, the 

velocity field (V) can be assumed small. Because of the 
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instability all the fields will be disturbed from their 
steady state values. If and 6P denote the 

fluctuations in temperature field, density and pressure 
field, then the onset of instability, they can be v/ritten 
as 


T = 

Ts 

+ 0 

(1.4a) 

P = 


+ 6P , 

(1>.4b) 

P = 


+ 6P • 

(1 .4c) 



where 



15 


a i Qj) 

m 

is the coefficient of thermal expansion. 

In deriving Eqs,(l .5-)-(t .7), Bossinesq approximation 
is applied. For variations in the temperature upto lO^K, 
the variations in the density are negligible because of 
the smallness of , The other coefficients (J'yA ) can be 
assumed constant. But the force of buoyancy due to variation 
in the density cannot be neglected. The magnitude of the 
buoyancy force is 


(6P) g = (p- Pjg = Pj^aCAT)g, (1 . 9 ) 

and this can be quite larger than other terms. Thus P 
(or Pg) is assum.ed a constant equal to P^^^ in all the terms 
of Navier-Stokes Eq,(ii,'7), except the term representing the 
force due to gravity. 


The linearized hydrodynamical ecyations 
for small disturbances are 


3 V _ 
3"T “ 


^6P 


- g a 0 z + V, 


m 


(t.10) 


“|| + (V.vTg) = 0, (1.11) 

'lAhere 

6P = _ aP^ 9 and g = -g2 (l .12) 

are used. 

The boundary conditions on temperature fluctuation 

are 
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0 at z 


= 1/2 . 


V At the rigid surfaces, not only the vertical component, 

but also the tangential components of the velocity field 
vanish because of ’no— slip’ condition. Using of this fact 
in Eq.(1 ^ ) yields 

If = 0 . (t.13) 

The term proportional to v5p in Eq,(1 ,7) is eliminated 
by taking curl of this equation twice. The z-component of 
the resulting equation is 


g 2 2 2 

V w = gav^G + J^Vw, 




vhere 


,2 = 
1 


3 
ax 


3^ 

— 2 
3y 


(1 .14) 


(l .15) 


To see the stability of the flovv^, the response of the system 
to all possible disturbances has to be investigated. This is 
done by expressing any arbitrary disturbance as a superposition 
of possible noimal modes and examining the stability of the 
system with respect of each of them. For the geometry 
considered here, the normal component (z-component) of 
velocity field and temperature fluctuation can be written 
in terms of two-dimensional periodic waves as given below: 

i(k^x + k„y) 


w = Vi?(z,t) e ^ ^ 

i( k X + k y) 

0 = 0 ( z, t) e ^ 


(tJ6) 
(t .17) 
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/ 2 2 

where k is the two-dimensional wave number 

associated mth the disturbances W(z,t} and 0(z,t), Using 
Eqs.(l.16), (1.1?) and (1. .S) in Eqs ,( ■} ,1 4) and(l..1l), 


9 , d^ 


k^) v! 


-g a +v 

dz 




W (^) + ^ * 

dz 


(1 .19) 


For convenience, the system of hydro dynamical equations 
is non-dimens ionali zed. The length is measured in the units 
of the separation 'd' between the plates, the time in units 

n 2, I 

of d /v , The non-dimensional wave number is defined as 
a Then Eqs.(l .1g) and (1^,19) become 

(D^-a^)(D^-a^ - Vf = 

and 

(D -a^ - p^ 0 = -(■^ d^) W , (tJ,2t) 

v;here, 

Pi (= “) is the Prandtl number 

'A 

and 

At the onset of stationary convection, the time dependence 
can be ignored in Eqs.(l ,20) and (1 ,2l). Thus, the relevant 
perturbation equations in non— dimensional form are 



8 


( 1 . 22 ) 


,^2 ^ 2.2 

( D - a } w = 


and 


/ —2 2\ '“v 

(D - a )(^m) = 


'•'7.1 th boundary conditions 

0 = W 


-Ra^ W 


= Di’ = 0 at z = 1/2. 


(1 .23) 


(1 .24) 


1 ,3 Chandrasekhar’s Technique 

The system of differential Eqs,(l.22) - (1 ,23) with 
boundary conditions (1,24) is solved using the following 
procedure. 

The temperature fluctuation field (0), on viiich the 

only requirement is the vanishing at the surfaces at 

z = + 1/2, is expanded in Fourier cosine series as 

oo ' 

0 = E \ Cos( 2n+1 ) 7 i;z. (1,25) 

n=0 

The above expansion is inserted in Eq,(l,22), This gives 

(D^-a^')^ VI = S cos( 2n+1 juz , (1,26) 

The general solution of this equation is 

cos( 2n+1 )%z 

= B cosh(az) + Cz sinh(az) + S „ , 

; {(2nt1 + a 1 

(1.27) 

Constants B and C are" determined from the constraints on 
W, They ‘are found to be 
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2B = -C tanh(a/2) 


(1 .28) 


and 


C = 


o 


2n+1 )7iA 


n 


(sinha-Pa J n ^ ^ ^2 ^2““2'j2' 


c-osh( a/ 2) , 


(1 . 29 ) 


Inserting Eqs.(l »25) and (1,27) in (l ,23) yields 


(D^-a^) ^ cos(2n+l)icz = 


— Ra rBcosh(az)+Cz sinh(az) +- 
cos( 2n+1 )7i;z 

+ ^ ?5”75 T’-!? — 

C(2n+1 rn^ + a^}'^ 


\'i^.'iich implies 

2 2 9 

{ ( 2n+1 ) Tt a ' 1 cos( 2n+1 )% z 

-l- 2 

n 


Ra^B cosh(az) +Cz sinh(az) 
Gos( 2n+1 ) z 
{(2n+l)^n^ + 

(1 ,39) 


Multiplying throughout vdth cos( 2n+1 )'n;z and integrating 
over the entire range of z( z = -1/2 to z = + 1/2), we 
find 


[C (2n+1 + a^i 


Ra' 


■ ' ^7~TT 

i{ 2n+1 ) % 




+ a-"} 


'] A. 


n 


+ 1/2 

= 2Ra^ B cosh(az) cos(2n+l)itz dz 

- 1/2 

+ 2R-'>^ C z sinh( az) ,cos( 2n+1 )'n;z dz . 

- 1/2 


( 1 ,3 1 ) 
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•Performing integrations and using Eqs, ( 1 , 28 ) and 
(1 . 29 ), 


[^( 2 n +1 ) } - 


Ka 


' ' ' ' ' ' " '} ') ^ 3 

{(2n+l)'"7T; + a } 


= - Ra [ 


16 a 


( 2 n-M )^‘n:Sij^dosh^( a/ 2 ) 


(sinha + a) 


C( 2 n +1 ) 




1 6 a oos h( a/ 2 ) ^ 2 m+ 1 ) ( 2 n +1 ) 


J 

(sihh a+a) '■ C( 2 m +1 2 n +1 )\^-+-a"^y 


(1.32) 


Now truncating the expansion (l., 25 ) at the first term 
( n=0) , 

,(o) 2 


,2 2,3 

+a ) 


[1 - 


16 7 ^^ a cosh^(a/ 2 ) 


( 'n;^+a‘^)^( sinha + a) 


], 


v^hich results in 


R 


(o) _ 


2 , 3 , 2 


( tt :^ + a }/a 


167 ^^ a cosh^(a/ 2 ) 


^ ^ r?+a^)^( sinh 


afa ) 


{^, 33 ) j 

' i 


To find critical value of R^ (= R°)> at v\hich the first 
convective instability occurs, R^ (a) in Eg. (l.. 33 ) is 
minimized with respect to ’a’. This yields a^ = 3,12 and 
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= 1712, within 1/j of the exact numerical results. 
Thus, the technique in the lowest order produce quite 
accurate results. The accuracy can be further improved by 
considering higher order terms in the expansion of temper- 
ature fluctuation field. 
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CHAPTER II 


CONVECTION IN BINARY LIQUIDS 


2.1 Introduction 

1 ~1 3 

The convection in a binary liquid — a mixture ' 
of two completely miscible and non-reactive liquids - — is 
more versatile than the single— componerat Rayleigh-Benard 
convection because of the possibility of both stationary 

and oscillatory instabilities at the onset. A mixture of 
3 4 

He- He near the superfluid transition temperature is 

the most general binary liquid mixture from the experimental 

point of view because of the follovang reasons; 

V 3 4 

i) the mass diffusion of He into He increases 

rapidly near the lamda point (superfluid 

transition temperature) and ultimately diverges 

.14 

at the transition point > while the themal 
conductivity at zero mass current remains 
finite; 

ii) the thermodiffusion -'coefficient, unlike that of 

any other binary liquid, becomes large and 

15 

reaches a limiting value of about 0,57 at the 
lamda point; and the effect of concentration 
gradient producing a heat current no longer be 
ignored; 
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iii) the coefficient of thermal expansion changes 
sign near the larada point ‘ , 

Because of the striking properties (i) and (ii),the 

ratio of mass diffusivity to the thermal conductivity (in 

3 4 

the absence of mass current) in the He- He mixture can be 
parametrically varied over a vd.de range — from a number much 
less than unity away from to a number much greater than 
unity very close to , Experimental advantages like good 

temperature stability and resolution provide a good control 
over the temperature field. Thus by varying the mean temper- 
ature of the liquid mixture all different possibilities at 

1 7—20 

the onset and subsequent transitions can be investigated. 

The onset of convection in binary liquids using rigid 

boundary conditions is discussed in depth by Gutkovdz— K3:usin 
21 

et al. They use a variational principle for stationary 

instability and wavefunctions similar to the variational ones 

for treating the oscillatory instability for which there is 

of 

no variational principle. In this chapter, the onset/convective 

instability — both stationary and oscillatory — vdth realistic 

boundary conditions has been studied using Chandrasekhar’s 
22 

technique , 

2,2 Hydrodynamics of Binary Liquids 

Hydrodynamic equations for the temperature (T) and 
concentration (®) variables in the binary mixture at the 
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onset have been obtained by Landau— Lifshitz as 


a T 
at 

ap 

ai 


(V,^) 

2 

= X V T , 

(2.1 ) 

(V.Jc) 

r 2 ^T 2J 

= Jd V C + — V 

(2.2) 


m 


because the heat current, due to concentration gradient, is 
negligibly small in most of the fluids. 

Here, V(X, t) is the velocity of an elementary volume 
of the liq-uid mixture at any point X in the system. The 
concentration ’c' is the mass fraction of He, oS", the 
isothermal mass diffusion coefficient, X , the thermal diffusi- 
vity in the absence of concentration gradient and k_, the 
thermodiffusion coefficient. 

In Boussinesq approximation, the Navier-Stokes equations 
describing velocity field in the liquid are given by 
-► 

^ + (V.v) V = - ^ VP + g -^ + a;v^V, (2,3) 

where is the mean density, P , the position dependent 

density and v , the kinematic viscosity. The incompressi- 
bility mf the fluid is expressed as 

V.v = 0 . (2.4) 

The equation of state is 

p> [i _ a(T - T^) - p(c-Cj^)3, (2.5) 

v\tiere Tj^^ and c^j^ are the mean temperature and concentration 
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respectively, a 


is the coefficient of thermal expansion 


- _ 1 3P 

a - • “ “pr“ -pf 

ra 


(2,6a) 


and 



"I dP 

~ T" aT' 

m 


(2.6b) 


measures the change in density of the mixture with change 
3 

in He concentration, F is positive according to the 
definition. 


In the steady conduction state, there is no motion 
in the fluid, i,e,, V = V = 0 and the temperatuie and 

o 

concentration fields are governed by the different equations 


= 0 (2.7a) 

3 z 


and 


3 2 c, 


s _ 


3 z' 


(2.7b) 


vliere the z— axis is positive in vertically upward direction. 

If and T 2 > ^2 temperatures - concentrations 

(of one component) of the mixture at the lower and upper 
plates respectively, steady state profiles for those variables 


are 


Vz) 

= T, - (A t) I , 

(2.8a) 

z) 

= - (A c) f , 

(2.8b) 

v#i ere AT = T| “T 2 

and Ac = - C 2 * 



As there is no mass current in the steady state, 



0 


Ac + ^- A T = 

that, is, ^ _ __T 

AT ~ • 

m 

If 6P , 6T, 6c and 6P be the fluctuations from steady state 
values for density, temperature, concentration and pressure 
because of the onset of instability, then hydrodynamic 
equation for these fluctuations are 


dV _ 
■ 


Pm 



( 2 . 10 ) 


+ W (2.12) 

m 

The teim v(6P) in Eq,(2,10) is eliminated by taking curl 
of the equation twice, as done in Chapter I, The z-coraponent 
of the resulting equation is 


w = g av^ (6T) + w , (2,13) 

2 3 ^ 8 ^ 

with V .. = — w + — 2 • 

“ 8 X 8y 

Now, Eqs ,( 2,1 1 )-( 2,13) are analyzed in normal modes by seeking 
solutions of the form 


6 T 

i( k X + k y) 

= @(z,t) .e , 

(2.14a) 

6 e 

i(k X + k y) 

s C ( z, t) . e ^ , 

(2,14b) 

w 

i(k„x + k y) 

= W( z, t) . e ^ , 

(2.14 c) 



where k ~ ^ wave number of the caivection 

rolls. 

By using Eqs,( 2,14a) - (2,14c) and Eqs .( 2.8a )-( 2,8b) in 
Eqs ,( 2,10 )-( 2,12) and linearizing the resulting equations, 
they can be written in non-dimensional form as 


(D^-a^)(D^ - - p) W = Ra^ (0- k C) , (2.15) 


( D^-a^ - ap )0 = -W , 

(D^-a^- |p)C= (D^ - a^)0, 


(2,16) 

(2.17) 


V . 

\Ahere p is the growth rate of fluctuations, cr( = -r) the 
thermal Prandtl number, R [ = — Rayleigh 

number, s the ratio of mass diffusivity (t>0) to thermal 

hf "g 

diffusivity (^), k(= ?f- §-) the separation parameter. In 

1 CY. 


'm 


Eqs ,( 2.1 5)-( 2,1 7) , the length is measured in units of 
separation (d) between the two plates, time, in units of 
d /i) , Q, in (AT), and C, in (Ac) . a =CH#d) 

Defining a new variable, the fluctuation in mass 

current 

Sj = (6c + 8T)=Aft [|f + ^) l|T)] = ic (C-gi), 

(2.18) 


m. 


and 


= = (c-gp. 


Now, from Eqs ,( 2,5) and (2.19), 


J 


(2.19) 





6P 

m 


= aA T [( 1-k)0 - kJ ] , 


( 2 . 20 ) 


Using Eq. (2,19) in Eqs ,( 2,.15)-( 2,17), the hydrodynamic 
•equations are transformed as follows; 

(D^-a^)(D^-a^-p)W = Ra^(l-k)0 - Ra^ kJ, 


(D^-a^- Cp) 0 = -W , 


(D^ - a^ - I p) J = (D^ 


2s © 

a ) - 

^ S 


( 2 , 21 ) 

( 2 . 22 ) 

{ 2 ^ 23 ) 


2.3 Boundary Conditions 

( i) Rigid bounding surfaces imply 
W = DW = 0 at z = + 1/2 


(2,24a) 


(ii) Since, the temperatures of lower and upper plates 


are fixed- and therefore, 

0 = 0 at z = + 1/2 


(2.241:) 


(iii) Impermeability prevents any mass current across 
the boundaries, that is, 


DJ = 0 at z = + 1/2 (2,24c) 

2,4 Stationary Convection 

In this case, the fluctuations are time— independent 
and, hence, p = 0, Now, following Chandrasekhar's procedure, 
the synmetric solution to the temperature fluctuation, 0, is 
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expanded in Fourier cosine series as 

oo 

0 = 2 A Cos(2n+l)nz, (2.25) 

n=0 

v\iiich is consistent vdth the boundary conditions on '0* . 
Inserting the expansion for temperature fluctuation from 
Eq,(2,25) in Eq.(2,23) mth the proper boundary conditions on 
mass current [Eq,( 2,24c) 1, fluctuation in mass current is 
determined exactly. It is found to be 

J = ^ [ n ^n GOs(2n+1 )'n;z + ^ Aj^( 2n+1 ) 

(2.26) 

2 2 

Operating by (D -a ) on Eg, (2,21) and using Eq.( 2,23) with 
p = 0, we find 

(D^-a^)\ = Ra^d-k- |)(D^-a^)©. (2,27) 


cosh( az ) 
sinh 




The additional conditions on the z-component of velocity 
field are 


(D^-a^)Sil = — kRa^J at z = + 1/2 , 


(2,28) 


Inserting Eg, (2, 25) in Eq,(2,27), the general solution of 
the latter is found to be 


2 

W( z) = c osh(az)+ B2Z sinh(az)+ B^z cosh(az) 
2 y ^n s^ 

+ Ra „ K — 55- TT-ry c •os( 2n+1 )it z 


(2.29) 



The boundary conditions in Eq,(2,28) yields 


8s 


Sinh ( a72) 


1 ^(-)"(2n.1)A„. 


(2.30) 


The boundary conditions W = DW = 0 at z = + 1/2 give 

B.] cosh(a/2) + sinh(a/2) = “ ~ Cosh(a/2) (2,31a) 


and 


a B.| sinh(a/2) + B 2 [sinh(a/ 2 ) + cosh(a/2)] 
= -BJ[^ + coth(a/2)3 einh(a/2) 


2 (-)’^(2n+l)A 

-r Ra (1— ir *^) 2 


(2,31b) 


B.| and B 2 are obtained frc«n the above two equations and 
are 


B 2 = + ^1 n^n *=oth(a/2)] dosh(a/2) , 

(sinh a+a) 

(2.32a) 

and 

®1 ~ ~ '5 ^ ®2 €^°'th(a/2) + “^] f (2,32b) 

v\tiere 


“1 = Ra^ (1 |), (2.33a) 

6 = ^ IL 

1 8s * a » 


(2.33b) 



n 


= C-)’^ ( 2n+1 )'n: A 


n. 


U2n-M )^ic^ + 


(2.33c) 


and 


n 


Qn = (“) (2n+l) A 


n 


(2.33d) 


Turning now to the differential Eq,( 2.22) [with p = 0 1 and 
using the solutions for and V*, it can be vnritten as 


2 2 

(D -a ) ^ cos( 2n+1 )nz = 


- B^cosh( az)-B 2 Z 6inh(az)- 


'B^z Cosh(az) - X- 


A 

^2" . ” "~2 yg" ^°s( 2n+1 )7Cz . 


C( 2n+1 )^ + a 

(2.34) 

Multiplying throughout with cos(2n+l)Ttz and integrating from 
2 — —1/2 to z = + 1/2, the above equation yields 


[(2n+l)^7t^ + a^] A = -^2 + 2B. x 

"" U2n+l)St2 + a^5^ ^ 


a. A 
1 n 


x^^^^cosh( az) cos( 2n+1 )'n:z dz + 

- 1/2 

+ 1/2 

+ 2Bo / 2 sinh(az), cos( 2n+1 )‘n :2 dz + 

-1/2 

+ 1/2 2 

+ 2Bo z c'osh( az) ,6os( 2n+1 )n;z dz 

- 1/2 

(2.35) 


Performing the integrations. 
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[!(2n+1)2n2 3^ j_ 1 ] 

{{2n+1)^n^+ 

4B^ ( -)'^ ( 2n-!-1 )7tCosh( a/2) 

{(2n+l)^ir^ + a'^3- " 


+ 


2 B 2 (-)’^ ( 2n + 1 ) TC 
{ ( 2n+1 + a^} 


x[s inh( a/2) 


4a cosh(a/2) 

— “ 

{(2n+1 ) Tc + a 


B3(-)’^(2n+1 ) 71 

2 - 

2{(2n+l) 7t + a } 


x{<iosh(a/2) - 


8 

{(2n+ 1)^11:''^ + 


2 

32 a cosh (a/2) 
{ t, 2n+1 ) Tc + a } 


Using Eqs.(2.30), ( 2.32a )-( 2.32b) and (2.33a)-(2 
above equation can be vnritten as 

^nn '^n m/=n ^mn '^m ~ ^ 

where 

^nn " ^ n C ( 1-^ |) V ^ ^ I 

^n 


- ]- 

} 


■{ C osh( a/2)+asinh( a/ ^ 


(2.36) 


.33d), the 


(2.37) 


(2n+l)^^j^(a), 


(2.38a) 
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1671^3 (2n+l)(2iiH-l)Y ^ cosh^( a/2) .Ra“( 1-k- -1) 


mn 


( Sinh a + a) 


+ |Ra^ (2n+1 )(2in+l) tt^Y ^ Dr.(a)» 


n 


^( 2n+1 + a^- } 


%(a) = 1 


1671^ a ( 2n+1 Cosh^C a/2) 

M„i ^ . -,. i_ . r« — nTl - i . 1 ii . . ■ .. »" ■ < ■ < , i . . ■■ .■i— 

( sinh a +a ) 


and 


a^Dn^a) = 4 Yj^ a^ Coth(a/2) -coth(a/2) ~ 


- a + 


2a coth( a/2) .cosh (a/2) 
(s'inh a + a) 


(2.3eb) 
( 2.38c) 

(2.38d) 


(2.38e) 


Eq, (2,37) is solvable if the determinant vanishes, which 
forces the value of R, the onset point for stationary convection. 
Thus R is the root of 


Det L^j = 0 


(2.39) 

If the expansion for ^ in Eq,(2,25) is truncated at the 
lowest order (i.e.,n = 0). 


Cos n: z 

and Eq,(2,37) gives 


Loo Ao = 0 • 


(2.40) 


From Eqs .( 2,38a )-( 2, 38e), and Eq.( 2.40), the value of R = (R^_°^) 





in the lowest order approximation is found to be 


>( o) _ 


-i- 


[ (1 -k- |)G^(a) + I G2(a).] 


(2.41) 


where 


and 


G,(a) = A ^(a) = 1 


1 6'n:^a cosh^C a/ 2) 

— — — , ■ - 

(■Jf'^+a ) (sinh a + a) 


( 2..42a) 


Gjla) = -u\{a) = [1 *4=^ 

ct n +a 


f + c osh a 
sinh a + a 


coth( a/ 2) ] 


(2..42b) 


The result inEq,( 2.40) is in agreement with one of the 
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variational approximation of Gutkowicz— Krusin et al. 


In most of the binary liquids and over a significant 

’ is very large compared 


3 4 

range in He— He mixture, 


to s. In this limit (k/s>>l), the value of is minimum 

s 

for a -► 0, The minimum value is given by, 

>( o) -2 2x3 


kR" 


~ tim 
^ 0 


+ a )' 


a^G^t a) - G^ a)3 


( 2.43) 


As, 


lim a^G.(a)= lija a^A(a) - 0 

a-^ 0 ' a 0 


( 2»43a) 
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- 2 ^2 ^ , _ 2 %^ - 24 


a G 2 (a) = - ti iin a'^D (a) 


a •*• 0 


Eq, (2,24) gives 


a -* 0 


(2.44b) 


a -♦ 0 


24 - 2 % 


^ 677 . 


(2.45) 


The exact numerical value for R is 720, Thus, the maximum 
error that one encounters in the lowest order approximation 
is only about 7%, The inclusion of the next order term 
(n=1) in the expansion of @ leads to the condition 


Loo h 1 


Lo1 ho - 0 . 


(2.46) 


In the limit k/s >> 1, according to Eqs.( 2,38a) - (2,38b) 
[as a O], 

I'ooHl - YoY? [(3 I p4’>)^f’-:VD„(a) D^(a)! 


k p( 1 ) . a^ ^ ^ ^ ^ L^^a) ^ 

s ^ 73 ^ 2 

Yo 


— -1 ] 

3 3-' 

Y oY 1 


(2.47a) 


Lolho - Y qY ^ Do^a) (2.47b) 
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because Xim a^A (a) = 0 = lim a^A.(a) . (2.47c) 

a 0 ° a 0 ‘ 

Substituting Eqs.( 2.47a) - (2.47b) in Eq. (2.46), the 

expression for ^ is found as 

s 


As 


S S 


0 = 


lin 
a 0 


1 

v2 ^3 r T^a^ 

0^1^ 3 + ' 

^ 1 o 


liia D.(a) 

a - 0 ° 



(2.48) 


Ha it^a^D/a) - , 

Eq, (2,48) yields 

- f -0 = 5 84^ ^92 -1 . - ? - - ^2,49) 

vhich is within ^% of the exact numerical answer. The results 
l-Eqs-V (2,4b) and (2,49)1 show the existence of stationary 
convection for zero wave number at the onset. That is, the 
wavelength is infinite, which means there is only one convection 
roll in the system at the onset of stationary convection, 

2,5 Oscillatory Convection 

In this case, the fluctuations are time dependent and 

Cii Q be the frequency of oscillation 


theyvary periodically. If 



at the 6nset then the velocity field and fluctuations in 
temperature and mass current can be written as 

id r 

W = W^gCz) e ® , (2.50a) 

© - ©os^^) ® ^ (2.50b) 

iS 

^ ® ° * ■ .^^*50g) 

So, for the case oscillatory instability, ’ p' is replaced by 

in Eqs .( 2,21 )*-( 2,23) and the relevant hydrodynamic equations 

ares 


(D^-a^)(D^-b^)W„ = Ros®^0-k)©„3 - - k J , 


= -«os* 

CD^-bg) = (D^-a^) 

v\tiere 

V.2 _ 2 

b — a + icd^ , 

2 2 — 
b^ = a + icrcd^ , 

I u2 2 .a — 

and b 2 = a + i — cd^ * 


© 

s * 


(2.51) 

(2.52) 

(2.53) 

( 2,54a) 

( 2,54b) 
( 2.54c) 


The boundary conditions remain the same as in Eqs, (2,24a) 
(2,24c), Since Chandrasekhar’s procedure yields quite good 
results for the threshold of stationary convection in the 
iowFest order we begin the study of oscillatory convection by 
making an ansatz on temperature fluctuation field as 
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e. 


os 


cosrtz , 
os 


(2.55) 


Inserting the above expression in Eg, (2,53) and thensolving 
it with boundary conditions DJ = 0 at z = + 1/2, we get the 
exact solution for the fluctuation in mass current given by 


os 


os 


(^^ + a^ ) r 
+ b|) 


c osh^ ^ 
^2 sinh(b2/2) 


+ COSTl z], 
(2.56) 


Now operating by (D -b^) on Eq,(2,5l) and using Eqs,(2.53) 
and (2,55), we have 


(D^-a^)(D^-b^)(D^-b|)v»^^ = 


R,3.A(1-k)(D2-b2)ei^-Ra%^ 

tD2-a2)0o^, 


-R a^[(1-k){n^+b2)'=ositz - I 


2 2 

("n; + a )] costiz 


■^os ^ ^os ’ 


(2.57^) 


v/here 


a 


= [(1-k)C®^ + bh - + a^)] 


[(1 -k- a^) + (1-k) • (2-5*) 


The differential Eq, ( 2.56) is to be solved with boundary 
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conditions given in Eq. (2.24a) and 


(D^-a^)(D^--b^-)W, 


P^os^ k J^3 at z = i 1/2 

(2.58) 


The general solution of Eq, (2,57) is 


C. cosh(az) + cosh (bz) + C,, cosh(b^z) 


R a a A ■ qos % z 
os os 


+ a^)(7^ ^+b^)('n; ^+b2) 


(2.59) 


The boundary conditions in Eq, (2,58) yields 


C = a 

3 os os 


2 k "n: 




(^ ^+b|)(b2-a^)(b|-b^)| 


(2,60) 


Other boundary conditions (W = DW = 0 at z = + 1/2) determine 
other constants, Vi/hich are found to be 


b^s inh( b2/2)cosh( b/2) - b sinh( b/2)cosh( b2/2)lC2 


“ V c»=''(b/2) _! 

( Tii'^+a )(^ +b )('^ "^^2 


(2.61a) 


>. ui /o^ cosh(b,^/2) 

, C o sh( a/ 2; , ^ 2 

'1 cosh(‘b72) ^3 cosh( b/2 ) •* ’ 


(2.61b) 


vi/here 



= [b Sinh( b/2)Cosh( 3 / 2 )- a SinhC a/2)Cosh( b/2) ] (?,61c) 


Now, inaerting Eqs,(2,55) and (2,59) in Eq,C2,52) and multi- 
plying throughout with cosstz in the resulting equation and 
then integrating from z = -1/2 to z = +1/2, we find 


{'^^+ b^) +1/2 +1/2 

= C, J" cosh(az) c-os'^z dz + c osh( bz).c os7t z dz 

2 ^ - 1/2 - 1/2 


+ 1/2 

+ Cq / cosh(b^z) oosTtz dz + 
- 1/2 ^ 


^ 2. CO 

.2 , 2\ /„ 2 , , 2\ /„ 2 


2(ix +a )(tc ^+b )(7i +bo) 


(11 ^+b^) 


■Alt [c 


.. cosh(a/2) cosh(b/2) cosh(b 2 / 2 ) 


1/ 2 2. 

( Tc + a ) 




V ' 1 

4. Cq . ---j- - - -j 

® + b?) 


os os 


(7t^ + a^)(7i;^ +‘b^) (it^ + bS) 


(2.62) 


Substituting f or , C 2 and from Eqs .(2,61a), (2,61b) and 

(2,60) respectively in the above equation and simplifying, we 


obtain 


(■re^+b^) 


A n -. 2 r 4 ^ k 
^os ^OS^ I- . 

s lo^ 


drah(b2/2) (i^^+bpCbl “ a ) 


2 9 9 

4 ^ (b -a )c,c>s_h 


(“rt^ + a^) 


— cosh(a/2) - — 2 COS h ( 62 / 2 )] 


(n + bo ) 


.cos( b/2) ^ 


(ii:^+a^)( Ab^)^(i^^ + b2) * (7c^+a^)('^^+b^)P^bi 


(2.63) 
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where ^2 - [ b 2 ®inh(b 2 / 2 ;oosh( b/2) - b dinh( b/2)cosh( b 2 / 2 ) 1 


(2.64) 


Aq ^ 0, the above equation yields 


^os 


+ b^)/a^ 


- ' -V n ^ 


[(1 >k- k/s)G^ (a,w) + IG^C a,w)] 


(2.65) 


where 


r 2 ,, 2 2 


{ -I _ .,1^. ■(,)? .-a .J g'^^hi .^’zL c^hi b/2). 


Gl(a-,9) - — 2 — 2 — "V — ^ - 

' ('n^^+a^)('rt -f-b^) A 


1 

(tc2+a^")(Tx2 


( 2.66a) 


and 


(a,«) 


471 ;^ 

BT- 


— 2 — o~^ — rj -- - 2 [ . A -I - Coth(b2/2) 

2 ( 7 t^+b|)(bJ - a^) ( 7 ^+b 2 ) 2 


- (b'-a") 


A 

“zrr 


2 c-osh( a /2) 1 (- ^ , ,, V 

^ sf^b^)"^ *■ (2.66b) 


In case k/s>>1 ^ the expression for oscillatory threshold 
becomes 


- = Re RT 

S OS 


/^2 2 n* 2 

+ b^)/a 


(2.67) 


[G^(a,a>) - a 2 (a, w) ] 

CN» 

where ^^(ajCo) is given by Eq,( 2,66b) with Vs >> 1, 

Thus, we find that the technique is equally applicable 
for both stationary as v;ell as oscillatory instabilities* 
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* 

CHAPTER III 


OVERSTABILITY IH K^GNETOHYDRODYNAMIC 
CONVECTIOM 


3.1 Introduction 

The effect of interaction between an external magnetic 
field and an electrically conducting fluid is two-fold; 

i) electric currents are generated because of the 

fluid motion in the magnetic field and, thus, the 
applied field is modified; 

ii) small fluid parcels carrying current experience 
an additional force (Lorentz force). 

1 

These effects, in general, delay the onset of thermal 

convection in the system. The possibility of overstability 

in the magnetohydrodynamic Benard convection is a problem 

of long standing. The problem was first analyzed by 
1 

Chandrasekhar , He used idealized boundary conditions to 
show that magnetic Prandtl number ( P 2 ) must be greater 
than thermal Prandtl number ( p^ ) for overstability to be- 
at all possible. He claimed that the principle of exchange 

of stability would be valid for P 2 < P-^ • Thus, he almost 

2 

ruled out the possibility of overstability. Gibson used 

realistic boundary conditions and obtained Chandrasekhar’s 

criterion for Q (Chandrasekhar’s number) >>1, Sherman 
3 

and Ostrach tackled the problan variationally and ,-f ot? the 
ease Q » 1 found the Chandrasekhar Gibson answer. Recently, 

* Contents of this chapter have been accepted in Phys . Fluids. 



4 

Banerjee et al, derived an exact condition, lAtiich states 

that the principle of exchange of stability is valid 

2 

provided Q p^ ^ % , This gives rise to an interesting 

possibility that for P2 < P-| j which is true for all practical 
cases, and for high enough magnetic field, it should be 
possible to obtain overstability. In this chapter, v^e 
investigate this possibility in RB geometry using realistic 
boundary conditions, 

3.2 Magnetohydrodynamic System 

A horiz0ni:ar layer of electrically conducting fluid 
confined between two rigid and perfectly conducting surfaces 
with a uniform magnetic field applied across it in vertically 
upward direction is heated from below. The equation of heat 
conduction remains the same as in the case of single -component 
fluid, but the equations of motion (Navier-Stokes equations) 
change due to "the presence of Loren tz. force. In Boussinesq 
approximation, they are 

|y + (V.v)v = - -pP + + 2 ^ V- V x H), 

Pm m 

(3.1) 

where 

J = o (E + X H) (3.2) 

is the current density in the fluid, a is the electrical 
conductivity and the magnetic pexme-abllityof the fluid, 

7he electric and magnetic fields are governed by Maxwell’s 



equations. For magnet chy dr odyn ami c systems, the displace- 
ment current is negligible in comparison to the induced 
current because of fluid motion, and Maxwell's equations 
are 


V • H = 0 , ( 3 ,3a ) 

V X H = J , (3.3b) 

V X E = — P-Qg ^ ^ (3,3c) 

The magnetic field is governed as 

-♦ 

^ ^ (V X H) = r\^ H ( 3 , 4 ) 

with 1/^^= ( crpo) • 

If h, 0 be the induced magnetic field, fluctuation in temper- 
ature field respectively at the onset of convection, then the 
z-component of Eq,(3,l) after elimination of vp from it and 
the equation for induced field are 


^ ^3t 


V ^)w = 


g a (V^Q 


329 

3^ 


and 


3h 


2 

^ h 




3 w 
3^ ^ 



^s 


(3,5) 


(3.6) 


where H is the magnetic field distribution in conduction 

3 

state. Now writing w and h as 


w = w ( z,t) 
h =h*(z,t) 


i( X + k^y) 
i(k^ X ‘+ k^y) 


f 


(3.7 a) 
(3.7b) 



“ ^“72 “ 2 “ 

where k = ’ is a two-dimensional wave number in 

XY plane. Now inserting Eos .( 3,7a)-( 3 ,7b) in Eqs,(3,5)- 
(3,6) and non-dim ensi on ali zing the resulting equations in 
linearized from, we obtain 


and 


(D^-a2)(D^-a^ - ^)^z,t) 

o o 3 -sf 

(D -a - P 2 “^)h (z,t) = 


= Ra^9 + Q(D^-a^)D h* ( z,t) 


Dv/ *( z,t) 


(3,8) 


where 

d^ 

Q = ~p — — is Chandrasekhar’s number, 

m o 

R, the Rayleigh number, d, the separation between the bounding 
surfaces and P 2 j the magnetic Prandtl number. 

Thus the relevant hydrodynamic equations for the 
simple magnetohydrodynamic Benard problem are 

(D^-a^)(D^-a^ - ^)w* = Ra^ +Q(D^-a^)DhT (3,9) 

( - a ~ P-] f ( 3 , 10 ) 

(D^ - a^ - P 2 ^)h = Dw . (3.11) 


3 ,3 Boundary Conditions 

The boundary surfaces considered here are assumed to 
be perfectly conducting. This means 

= 0 at z =+ 1/2 . 

Other boundary conditions because of rigidity of the bounding 
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surfaces are 

w* = Dv*= 9 = 0 at z = + 1/2 . 

3 .4 Mathematical Analysis. 

Since the analysis is intended to look for any 
possibility of overstability, the variables are expressed 
as 

/ w * \ == / ’'■?( z ) \ 

> ! \ icom 

I Biz) S e ° , (3.12) 

. K(z) / 

\ J 

where is the frequency of the convection roll at the 
onset of overstability. Therefore, Eqs .( 3,9)“( 3 ,1 1 ) become 


(D^-a^)(D^-b^)W = Ra^0 + Q(D^-a^)DK, (3.13) 

(D^ - e = -V*/, (3.14) 

(D^-Y^) K = DW, (3.15) 

where 

2 2 2 2 2 2 
b = a + ic<)o, P = a + and Y = a + im^p^. 

(3.16) 

The boundary conditions are now 

W = DW =@ = K = 0 at z = + 1/2 . (3.17) 



Following the procedure discussed earlier, the temperature 
field is expanded in Fourier cosine series as 



DU 


e = 

2 A c os ( 2n 1 ) nz , 

n=0 

• 

(3.18) 

Now-. E5(.it-3r;i4 ) 

is inserted in Eq,(3,l5), 

Vvhich is then solved 

using Eq,(3,1 8) 

and boundary conditions 

on K , The solution 

is given by 




K= 2 (2n+l)A [sin( 2n+1 )Ti:z 

n (2n+l)"^ + Y 


Inserting the expansion for0 and K in Eq.CS.Ia), the latter 
becomes 


[ (D^-a^)(D^-b^) - Q D^]W = [Ra^ + i Q(2n+l)^'rt^ 


X . 't . P . ,. — ] cos( 2n+1 )-jt;z - 

( ?n+1 ) % + Y 

-mp^TQ .A„(2n.i)n-)" 

V cosKy z) 

cosw/:^ • 


(3.20) 


The general solution of this equation is 

P. 


W = 2 A, 


1n 


n ^ { ( 2n+1 ) ^ + q 2n-t t )“ A q^ > U 2n+ 1 ) ^^+Y^ 5 


C(Ss( 2n+1 )7ii 


+ 2 Aj^ P 2 ^ cosh(Yz)4- B. cosh(q^z)+ coshCq^z), (3,2l) 

n 

vjhere 

P^ =[Ra^{( 2n+1 + Y^ > -H«oP2 ^ 2n+1 2n+1 ) 

(3.20 
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i«P2YQ(2n-i-l)jt(-)'^C(2n+1 ] 

{(2n+1)\^ -;-T'^}{qY^ - (Y^-a^)(Y^-b^) sinh(Y/2) 

(3.23) 

2 [iQ + (a+b)^5^/^ + £Q + (a-b)^}^/^], {3,2^) 

The boundary conditions on W give 

cosh(q^/2) + Cosh(q2/2) = - S P ^ cosh(|r/2) 

n ^ 

(3.25) 

and 

q^B^ sinh(q^/2) + q^B^ 3inh(q2/2) = - 2 Y sinh(]r/2 ) 

+ S P.n (2n+l) (-*)^ie (3.26) 

n 

leading to 


2n 


and 


h ,2 


■q^ sinh(-^) ^ A^ P^^ cosh(Y/2)+ \ A^ P2r^ Y sinHh/2) eosh( 


B. 


- oosh(q2/2) n \ P-|n(2n+1>(-)'^^] 




^1- 


[ q2 sinh ( ”) cosh ( q-j einh ( -“) cosh ( -^ ) ] 


( 3 . 27 ) 


and 


10B284 


B 2 = -[ 2 A^ P 2 j^ c osh(Y/2)+ B^ cosKc|^/2)l/ cosh(q2/2) 

■ (3,28) 

mth 

^ = P-| ^CK 2n+1 + q^H( 2rd-1 + q^ H( 2n+1 + Y ^ ^3. 

( 3 . 29 ) 
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Inserting the solution given in Eq,(3,2l) together with the 
expansion of 00,(3,18'), in Eg, (3, 14), multiplying both sides 
with (^&${2n+^)%z and integrating from z === -1/2 to z = +1/2, 
the following expression is obtained 

( ( 2n+1 }{( 2m+1 )^7t^ + } {( 2n+1 )V + p ^} {( 2m+l)«^+ Y 


A [ 1 - A ] - 

m 1 m *- 


mm‘ 


A P, 


A 


n/ro n In mn 


+ i \ P2„ L„[(2m+l)V 


(3,30) 


where 


mn 


4‘n^( 2n+1 )( 2m+1 o^-c^) cosh(^) cosK ) 
C ( 2m+1 )^TE^+q| H ( 2m+1 )^‘Jt^+q 2 einh( -^) GOsh{-^) 


q^ sint( -or) sinh("?f ) ] 


- 


(3.3f) 


and 


Ljn ^ ( •")’” 4C28H-1) ^ 


(q^ -Y ^)?(2n+l)^^ + q?> o>osH!(/2) 


C (2m+V)^it^ + Y^ 3 


4- 


4 


2 2 Y Y ^2 

( q|-q2) cos h(-5) C sint^ cosh(-^)~ q^ cosh(2)sinh(-^) J 

^ ^ - -q ' _ 

? q^ S'inh(^) Posh ( q2 S'inh(-^) cosh ( ) * 


( 3^32) 


Befining, 
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i ( 2m+1 -I- q^}£(2iiH-1 + q|}{(2m+1 + p ^ } 

i ( 2m+1 + Y p ( 1- A )_p l £(2m+1 ^ + Y ^ ^ 

1m ^ mm^ 2m m '• ’ 


and 


2„2 


ran 


= - P2n U2tn+1) it‘ + Y^!, 


(3.3 3) 

(3.34) 


Eq. (3,30) can be rewritten as 

Gmr, = 0 . (3.35) 

min m npm mn n 

The consistency condition is 

Det G = 0, (3.36) 

The entries of the determinant being complex numbers, yield 
two equations which are solved for R and 

In principle this yields the exact solution. In 
practice, the evaluation of large determinants is not 
necessary. As shown in Chapters I andil, the lowest order 
truncation yields quite accurate results. In fact one-mode 
(n = 0) truncation is extremely good for a single-component 
fluid. 

3,5 Overstability 

For overstability one-mode truncation should yield 
accurate results in the low frequency regime. At high 
frequencies the truncation is susceptible to error as larger 



determinants are going to produce higher powers of m^and 

lead to inconsistency. Now the analysis can be simplified 

2 

by considering the limit Q >>1, P 2 < < 1 'with Qp^C > ) 
finite. Equating real and imaginary parts of Eq4(3*36) in 
its lowest order approximation, one obtains 


(n: ^+a^)p^ P2 ^ (H-P^+P2 ) + Q ^+a^) x 

X ( P^ ‘''P 2 ^ “ ^2^ ”* ^2 ^ ^ •3-7} 

and 

Ra^ = + Q7i:^(Aa^) - P2^ -“I ) 

(it +a ) 

+ (ic^+a^)^ (p-j+P2^ P-|P2^)^» (3,3‘8) 


where 


A 


A- 


. 4a Coth( a/2) 

. 2 2,^ 
(it + a ) 




(3,39) 


For idealized boundary conditions A = 1 and the equations 
become similar to those obtained . by Chandrasekhar , It is 
the use of proper boundary conditions , vAich make ^ 1 


and lead to the possibility of overstability even though 
P 2 <<^1 • Eliminating R from Eq,(3,37) by using Eq,(3,38), and 
working to the lowest ordei in P 2 yields 



Po Q Pi + a^)^(1+P^) 

^ r + i-3, (340) 

P ’1 ( A - 1 ) (A - 1 ) 


where Q = 
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Since the wave number 'a' increases with 'Q’ and L - 1 
for Q >> 1 , the Q dependence of various quantities are 
self consistent if, in accordance with the case for 
stationary instability, a m In this case, the first 

term of Eq, (3,4 0) dominates and 

r \ 

Q Q P2 7 i; 2 /(a -1 ) p^it (a/4) , (3.41) 


It is interesting to note that Qp2 enters as a combination 
in the expression for the frequency of overstability. 

The Rayleigh number at this order of calculation is 


Ra^ 


^ (ir^+a^)^ + Q 71^(71^3^) - 


2 ,„ 2 , 2 ^ 


( 




Q7^P2)" 




■)• 


(3.42) 


which is lower than the threshold R^ for stationary 
convection given by 

2 9 9 q 9 9 0 

R5 a = + a ) + +a ). (3.43) 


Thus the possibility of overstability even if P2 < < P-j is 
not ruled out. 


3 ,6 Necessary Condition for Existence of a Polycritical Point 
If there is an overstable convection, it is relevant to 
ask whether there is a point in paraneter space where the 
lines of oscillatory and stationary convections meet. This 
will indicate the possibility of any polycritical point where 
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both types of phenomena are possible in the system. The 

condition for this possibility is obtained by equating the 

2 

expression for to zero. With one— mode truncation and 
in the limit Q >>1, this condition [from Eqsiip.S?)— ( 3,39) ] 
becomes 

Q Tc ( IX +a )(p^-P 2 ) + ( 1+p-j )( + 4nP2 ^ ® coth( a/2) 

= 0 (3.44) 

Since the last two terms on the l.h.s, are positive definite, 
the condition is fulfilled only if 

< P2 » 

which is a necessary condition for the existence of 

5-6 

codimension-2 bifurcations . Thus Chandrasekhar-Gibson 
criterion is seen to be a necessary condition for the 
existence of a polycritical point and not a necessary condition 
for overstability. 
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CHAPTER IV * 

CONVECTION IN MODULATED -RAYLEIGH-B&ARD FLOW 
4,1 Introduction 

In this chapter the onset of convective instability 
in a modulated Rayleigh-Benard system mth realistic 
boundaries is discussed. An external modulation of the 
temperature difference betv^en the plates can drastically 
modify the behaviour of the system. The study of the effect 
of modulation is currently receiving considerable attention, 
spurred by highly sophisticated experimental techniques, A 
proper tuning of amplitude and frequency of modulation can 
lead to a new class of phenomena and thus makes the problem 
of hydrodynamic instability more versatile and interesting, 

1—10 

The theoretical treatnent of the effect of 

1 —4 

modulation is based either on full hydrodynamic equations 

5-1 0 

or a few-mode truncation thereof , Vvhile the few-mode 
truncation always yields solutions in closed form, the 
hydrodynamic equations seem to yield closed-form solutions 
only v^?ith idealized (free) boundary conditions. Here, 
using Chandrasekhar’s procedure, the critical Raylei^ number 
is obtained in closed form for the system with proper rigid 
boundary conditions, 

* Contents of this chapter have been accepted in Phys ,Rev, A, 
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*2 jiZdrodYnamics of Modulated Rayleigh-Benard Flow 

The control parameter — ” the temperature difference 
between the two plates ““ is modulated by imposing an 
external time-periodic disturbance on the temperature of 
the lower plate. The temperature (Tj_) of the lower plate 
is 

T^ = T^ + eRe(AT) e^“^, (4.1) 


where e and m are the amplitude and frequency of the modula- 
tion respectively, T^ and T 2 are the temperatures of lower 
and upper plates respectively in the absence of any modulation, 
At = T^ - T 2 » Since the equation governing the conduction 
of heat is linear and is given by 


!!h 

at 



(4.2) 


the temperature field (Tjj) of the hydrodynamic system can 
be expressed as 


Tpj(Z,t) = Tg(Z) + eT^g(Z,t) . (4.3) 


'T ’ is the steady state temperature field in the absence 
$ 

of any modulations 3nd is given as 

^^(Z) = T^ - (AT)(^ + ^) . ' (4.4) 

If *Z’ is measured in terms of the separation ’d’ between 
the plates, the above equation takes the form 

T (z) = T. - (AT)(z + I ), 

s * 


(4.5a) 
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Th - - AT(z + I) eRe(AT) 


sinh {odC-^ - z)} 
sinh(a d) 


icot 

e 


(4.10) 


If ©, as defined earlier (Chapter I), be the fluctuation in 
the temperature field at the onset of convection, the 
modified temperature field will be 

T = + 0 . (4,11) 

Now the linearized perturbation equations, as obtained in 
Chapter I, are 

(D^-a^)(D^ - - Ip) w = 0 (4.12) 

and 

(D^-a^ ~Pi'It) ^ ~ + Re f(z) 

(4.13) 

where 

' si ' n^^d)" cosh{ad(^- z) }, (4.14) 

- ^ 

and Q, m and -x are measured in units of AT, ~ ~o '~ snd — 

d^ y 

respectively, p^ and R have their usual meaning^ As the 
realistic bounding surfaces are rigid, the boundary conditions 
are 

w = Dw = 0= Oatz=+ , (4,15) 


4,3 Mathematical Analysis 

The critical value R can be obtained by solving the 
Eqs ,( 4.1 2)-( 4,13) with boundary conditions given by 
Eq,(4,15), For small modulation ( e <<1), the solutions are 


t p is thermal Prandtle number and R, the Rayleigh number 
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obtained by expanding w, G and R in powers of ‘e ’ . That is, 



% 


e 

- vv^ 

4- 

' ♦ • af f 

(4.16) 

0 

= ®o 

4- 

e 

+ ©2 

♦ f 

(4,17) 

R 

= ^0 

4" 

e R.| 

+ R2 


(4*18) 


Inserting these expansions in Eqs ,( 4,1 2}-( 4,13) and (4,15) 

and equating terms of the same order in 'e' , the resulting 

2 

equations up to 0(e - ) are 


(D^-atCD^-a^ - - e,. 

(4.19a) 


(4.19b) 

(D^-at (D^ - -3T)®1 = ®1 ’ 

(4.20a) 

(D^-a^ -p, 4-^)9^ = - - 

„ 2 _ , _ iSr,, 

R^a Re( fw^e ) ; 

(4.20b) 

9 9 9 9 3 

(D -a )iD^ - a - ^)w 2 = ©2, 

(4.21a) 


(D^-a^~p^-|^e2 = -V% "■ ReCfw^e'' “^) 

-R^a^ Re(fwQ “^ ) (4.21b) 

The zeroth order equations represent the system in the absence 
of any modulation. In this situation the principle of 
exchange of stabilities is valid and these equations are 
reduced to 
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(D^ - 




w = 


(4.22a) 


and 



2 -^ 
a )Q^ 




(4.22b) 


The realistic boundary conditions are 

©Q = = DWq » 0 at z=+1/2, (4,23) 

The solutions will be symmetric at the lowest Rayleigh 

1 1 

number, Chandrasekhar’s technique is applied to solve 
these equations. As one— mode truncation yields extremely 
accurate results (Chapter I) for a single-component fluid, 
the temperature fluctuation *0' is expanded in Fourier 
cosine series and truncated at the first term. Therefore, 

cosiCz , (4,24) 

Inserting the above in Eq,( 4,22a), the solution for w^ 
with boundary conditions on [Eg. (4,23) 3 becomes 

Aq cos-rcz 

° ° ° (71 2 + 

(4,25) 

where 

2Bq = “• EIq tanh (a/2) (4,26) 

and / , N 

2tcA Coshva/2) 

C, = - 7% -^ — . (4.27) 

(Tt +a )'^(sinh +a) 

Inserting Eqs,( 4,24) and (4.25) in Eq,( 4.22b) , 
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multiplying throughout with cosnz and integrating from 
z = -1/2 to z = 1/2 yields 

R = ( r 1 _ 16^. cosh^C a /a), _ 

(4.28) 

which is the same result as obtained in Chapter l[Eq,( t ,J3)} 

The perturbative corrections etc,, vhich 

are compatible with may be obtained by requiring that 
Eqs ,( 4.20a)-( 4,20b) and Eqs ,( 4,21a)-(4,21b) must have 
solutions compatible with the solutions of Eqs ,( 2.19a)— 
(2,19b)« The solvability criterion for Eqs ,( 4,21 a)-( 4,21b) 
demands that 

0 \ 


1 , ^ rs>. , 

T ^ ®o^ 


dT^ = 0, 

\ -R a^w -R^a^Re( fw^e^ “'^)/ 

\ o o 0 ^ o } 


where w^ and 0^ ore left reactors of Eqs ,( 2,1 9a )-( 2,19b) , 


It leads to 


= 0 


(4,30) 


To obtain the forms of and 9^ within the one mode Fourier 
approximation described above, fw^ is expanded in Fourier 
series and only the costiz term is retained. We write 


f Wq ^ 


vhere 


F cosJtz , 


( 4,31 ) 
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+ 1/2 

F = 2 / f v\? c'osTiz dz 

- 1/2 ° 




'1'^ 


( “ { 4%\{a d )^ 


+ 


B ad 
sirf^d )' 


,71 Cosh C-^) 2 : 


x[ 


Coshli (a-ad)} coshCi (a+ad)l C 

I ^ T I -O 

— ■ fs M “T — " r * " ' : " """ " ; '?'’ " 


{( a-ad) + 11 } { ( a+ad)^+ 3 


2 “ r 2 T ^ ''’ sinhta d )^ 


Q,, sinh {x( a-ad)3 sinh{i ( a+ad)} 
XTt cosh(-7r) [-— -Q — r?5- + o 


C ( a-ad) + Tc 


{(a-ttxd) + % 3 


4( a-ad) 
{( a - ad )''^+ 


(a-ad) 4( a +ad) 

Cosh C ^rs ' ' 3 " Q o o ^ 

^ {(a +ad)^+it 3^ 




(4.32) 


iUT 


The time dependence of w.| and 8 .j is of the-form e , so that 
and G.J can be '.'written as 


®1 = 

^ f . iCi>T 

0 -I (z) e , 

(4.33a) 

W.| = 

Hyvz) e 

(4,33b) 

Substituting above 

in Eqs ,( 4.20a)-( 4,20b) with R.| 

= 0 , we 


find 

( D ^- a ^)( D ^- b ^) 
( D ^- b 2 ) e ^ = - W , 

wtiere, 

2 2 . - 

b = a + i u . 

2 2 — 
b 2 = a + io ) P . J . 


(0 ^ i ( 4,34a) 

- R^a^ fw^ (4.34b) 

(4,35a) 


and 


(4,35 b) 
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The boundary conditions are 

©1 = = DW^ = 0 at z = + 1/2 . (4,36) 

Bf pfoeeeding, as for zeroth order solution, by making the 
ansatz 

0^ = cosuz , (4,37) 

the solution for is obtained from Eq,( 4,34a) v.dth boundary 
conditions on in Eq,(4,36), The solution is found to be 


A. bos-iiz 

cosh(az) + C. cosh(bz) + — g — ^ — 2 — * 

(ti +a )('rc +b ) 

(4,38) 

vdth 


c’osh(a/2) = -C^ cosh(b/2) 


(4.39) 


and 

C 


cosh ( a/2) 

^ ^ [h sinh(b/2)c'osh( a/2)--a dinh( a/2)cosh( b/23 


1 

' ' ■ W - "VT" 
(n^+a^) 


1 


T 


(% + b") 

The amplitude A^ is obtained from Eq,( 4.20b) by using the 
above equation [ Eq.( 4,37)3 and Eqs ,( 4.38)-( 4,40) as 


(4.40) 


-R-a-Sr 


R F 
o 


(It ^+a^-fpdu) - ^ 

^ (ii^+a^)(ic^+b 


i' 


4^a^~lcr} cosh( a/2)cosh( b/ 2 ) 

3inh(b/2^ ; 

<^osh( a/2)--a sinh(a/2)- 
-fcosh(i)/2)} 


(4.41) 
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2 

Now turning to 0(e ) equations and applying solva- 
bility condition, we find that 

T — 

Y ^ [Rg® < ©Q ! vv^>+ <0^ ! fw^ e-" >]dT= 0 . 

(4.42) 

Inserting the expressions for w^ and taking the time average, 
we arrive at 


Ro = 


-Ro Re 


<0 1 f W. > 

. O 1 


= -R„ Re 


2c } w^> 
o * o 

< ^+a")©Qf > 

2< ( T?'+a^)Q^\ v;^ > 


= -R. Re 


<(D^-a^)0olf > 

2<( D^-a^) ©^ I w^> 


- R^ Re 
o 


<w I f W. > 
o 1 

2 <-'o • '"o> 


= - Ro Re 


< 'iL 1 f w > 

T O 


+ 1/2 * 

Evaluating <W^|fw^> = / f w^ dz 


and 


- 1/2 
+ 1/2 ^ 

< w.|w > = / w w dz 

o o ° 


(4.43) 

(4.44a) 

(4.44b) 


and inserting in Eg, (4.43), we find 


rI IFI^ 

R2 = “2— — Re 

4<WqIw^> 


G( G)) 


^ r ^ - , , , ^ 7 

RQa^G(m )-( % +a ) (it ^+a +i6) )( ic^+a + i P-, © ) 

(4.45) 
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where 


G(g)) = 1 + 


4it^( - b^) cosh( a/ 2) cosh( b/ 2) 

“^ 2 ^ 2 — '*"'2 — *" 2 " ' — ^ — ■ — - — — ■ 

(u +a )(-n; +b )[ b s inh( b/2) cosh(a/2)~a sinh( a/2)cosh( b/2)] 


(4.46) 


4 ,4 Results and Discussions 

It is relevant to look at the asymptotic behaviours of 
^2'^^o ' ~ fractional correction to the critical Rayleigh 

number because of external modulation. In zero frequency 
limit, G(m ) is expanded in powers of iwas 


G(m) = G(o) + (im) G'(S)j + (io>)^G"(w)| + 


• • • • » 


m=o 


0=0 


(4.47) 


where 


16it^a cosh^(a/2) 

G ( O ) = 1 ~ 9 9 9 ^ 

(^ +a ) (sinh a+a) 


(4,48a) 


and 


G(0)! « 


( + a^)^(Pl+ 1) 

^o ® 


(4.48b) 


and the second derivative contributes even less. 

In this limit, the fractional correction turns out to 
be 


-R— ^ - 2 [— » 


(F^ D' 


(4.49) 


1 


which is same as that with idealized boundary conditions . 
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That is for very small frequencies, the modulation stabilizes 
the conduction stats and correction to the critical Rayleigh 
number is independent of the types of bounding surfaces 
considered. 

In very high frequency ( m -► oo) regime. 


G(m) - 1, 


and thus 



I F( to) I 

a 





(4.50) 


The correction is positive definite and monotonic. It is 
different than the results for idealized boundary conditions, 
where the fractional correction is proportional to 5 ” , 

Now we compare our results with the relevant numerical 

1 

work by Rosenblat and Tanaka * Choosing the am.piitude of 
modulation e= 0,4, our results agree with that of the above 
mentioned numerical work for both p^= 10 and p^= 1 within 

b% for all frequencies. This establishes the soundness of 
the procedure applied in the present work. 
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CHAPTER V 

CONVECTION IN MODUL/vTED BIM^RY LIQUID S 

5.1 Introduction 

1-3 

RB convection in binary liquids has of late attra- 
cted a good deal of attention. The possibility of oscilla— 

4 5 

tory convection’ at the onset of hydrodynamic instability 
and subsequent codimensior>-two bifurcations has prompted 
the recent spate of experiments , Because of the above two 
possibilities, parametric m.odulation of the temperature 
difference between the two plates is expected to produce 
interesting effects on the flo’W pattern. In this chapter, 
the effect of modulation on the onset of stationary as well 
as oscillatory convection is investigated .for binary liquid 
mixture confined to realistic bounding surfaces, 

n 

Chandrasekhar's procedure is applied to incorporate the 
realistic boundary conditions, 

5 .2 Hydrodynamics of Modulated Binary Liquids 

Equations describing the fluid motion, temperature 
and concentration fields remain similar to those given in 
Chapter ii[Eqs .( 2.1 )-( 2.3) ] . The incompressibility condition 
and the equation of state are also given by Eqs .( 2.4)-( 2.6) . 

As the tempera tvire of the lower plate is modulated as 

Tj^ “ ^1 ^ (5.1) 


* Contents of this cliapter have been submitted for publication. 



the temperature and concentration profiles in steady state 
become. 


T3(z) 


T-j " (at) ('^ + 'i) + e Re{ 


(AT) sinh?i((-^ - z) }e^^ 


'ainh (yd ) 


(5.2) 


03(2) = 


o, + (AT) t; ® rriara) 




-CoshJT](-^ + z)}- coshClfd)coshC'n(^ - z)} 

©inhCnd) 


) - 0] 


(5,3a) 


where 


Y 



i^ 

K ’ 

im 

J3 ^ 


(5.3b) 

(5.3c) 


and other symbols have their usual meaning as defined in 
Chapter II, At the onset of convection, fluid motion starts 
in the ’system. If V re pP-T* S 3 nts velocity at the onset and 
6P ,6 T , 6c and 6p, the fluctuations in mean density, 
temperature, concentration and pressure respectively because 
of fluid motion, the hydrodynamic equations for disturbances 
in linearized form are given by Eqs .( 2,1 0)-( 2,12) in 
Chapter H, The infinite extensions of bounding surfaces 
allow us to express these disturbances in terms of two- 
dimensional periodic waves in XY plane. The hydrodynamic 
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equations on non— diraensionalization are expressed as 


(D^-a^)(D^-a2 - |^)w = Ra^ (©- kc), (^.4) 

(D^-a^- 0|-)0 =-w{l+ f(z, 6 ))e^“^}, (5.5) 

= (D^-a^)© -f[l - g(z,6>)e^“'"3,(5.6) 


v^h ere 




and 


g( z 


»co 


(5.7a) 


.r—-^ Tid sinl^. (z + 4)l( 1 +dosWc5) J 

Re. F— ^ 1 

(5.7b) 


By defining a non-dimensional mass current 


j = c - 9 (5.8) 

and dropping the correction in the differential equation for 
j, the linearized perturbation equations for the system of 
modulated binary liquids are 

(D^-a^)(D^-a^-p)w = Ra^(l-k)e - Ra^k j, (5.9) 

o o I 1 

(D -a^ - qp)© = -w[ 1 + f(z,w) e 4 (5,10) 

(D^-a^- |p )3 - (D^-a^)| , (5.11) 

where 'a' is non-dimensional wave number in horizontal plane 
and * p’ represents the rate of growth of fluctuations. 
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The realistic boundary conditions are 

w = Dw = 9 = Dj = 0 at z = + 1/2 . (5.12) 

We proceed by assuming the amplitude of modulation to 
be small ( e «1 ) , in this case, we expand w, 0, j and R in 


powers of 'e’ as 

w = + ew^ + w^ + ... , (5.13) 

0 = 0^ + 60^ + e- 9^ + ... , (5.14) 

j = in jg ••• * (5,15) 

R = Rq + eR^ + R2+ ... , (5.16) 


Inserting these expansions in Eqs .( 5 ,9)-( 5,1 2) and collecting 

2 

terms of the same order in ’e', we find up to 0(e ) 


( D^-a^ 

)(D^-a%)w^ 



ic)e„ - 

kJo> 

(5.17a) 

(D^-a^ 

- c7p)e^ 

=: 




(5.17b) 

(D^-a^ 

- 1 P)jo 

= 

(D^-a^) 

fo 

S 


(5.17c) 

with 







% = 

Dw = 6^ = 
0 0 


= 0 at 

Z = ± 

1 . 

2 ’ 

(5.17d) 

(D^-a 

^)(D^-a^-p)w^ 

= 

V^d- 

k)©-j - 

R^a^k 





+R^a^(1- 

k)9o - 

R,a2 kj„. 

(5.1 an) 

(D^-a 

dp)0^ 


-w^ - f(z,(i)e 

iirc 

Wo, 

(5.18b) 
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(D^-a^ - I p)j^ ^ (5.18c) 

wi th 


W.J = 

Dw.j = 0 ^ 

- Dj.| = 0 at z = + 5 

(5.18d) 

(D^-a^)(D^-a^- 

■p)w 2 = RQa"(l-k)e 2 - 




+ R.ja^(l-k)©., + R2a^(l-k)eQ 




- R^a^ kj., - R 2 a^ kj^. 

(5.19a) 

1 

CM 

cd 

1 

CM 

Q 

p )02 

-/ N iSc' 

= -w^ - f( z,u )e W.J , 

(5.19b) 

/ 2 2 

( D -a^- 

C'’ \ 

- 3 P) 

= (D 2 -a 2 )i 2 

(5.19c) 


vdth 

W 2 = Dw^ = 02 ~ Dj2 = 0a'tz = +-2. (5.19d) 

Equations ( 5,f7a) — (5J7c) represent unmodulated system 
of the binary liquid, mixtures. 


5.3 Stationary Convection and Modulation 

At the onset of stationary convection the fluctuations 
are time -independent ( p=0) in the absence of external modu- 
lation. So, Eqs .( 5,17a)-( 5.17c) reduce to 


(D^-a )\ = R^a^(1-k)6^ - 





f 


( 5*20a) 
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(D^ - a^)0^ = ^ (5,20b) 

(D^ - a^)jo = (D^ - a^-) (5.20c) 

with boundary conditions given by Eq.( 5 .17u) , To solve this 

system of differential equations, we follov/ Chandrasekhar’ s 
7 

technique and expand ©^in a Fourier cosine series as 


00 

©_ = ^ K ^os(2n+l)iiz , (5.21) 

° n=0 

As we have seen in Chapter II, one-mode truncation yields 
quite accurate results (within 7%) and, therefore, we 
determine the effect of modulation within one-mode approxi- 
mation. We choose - 


9^ = Aq costcz , 


(5.22) 


With this ansatz on ’9^, the Eq.,(5,20c) 
boundary conditions Dj^ =0atz=+1/2 
is found to be 

_ ^ cosh( az) 

•^o s " a ”s'inh( a/2; 

2 . ' 

Now operating by ( D - a ; on Eq,( 5,20a) 
the former becomes 


is solved with 
The solution 

• cositz], (5.23) 

and using Eq.,(5..20c) 


f 


(D— a^)\ = (1 -k- |)(D^-a^)0^ . (5.24) 

As 0 - 0 at z = + 1/2, Eq,( 5.20a) yields a boundary condition 

O 

on Wq given by. . 
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2^2 O ^ 

® ) % 7 “VIOq at z = + 1 (5,25) 

in addition to the boundary conditions v/ = Dw = O' at 

o o 

z = + 1/'2. Inserting Eq.(5.22) in Eq.(5.24), the latter 
bee omes 

2 7) 3 o i_ 

(D -a^) w = - P a (1 _ ic- + a^)A costez , (5.26) 

U o o 

The general solution of this equation is 


w. 


“ ^o ‘-osh(az) -i- B 2 Z ginh(az) + cosh(az) 


+ Rq a 


(1 - k- i) 

Q 

( ^^ + a^)^ 


COStiZ] , 


(5.27) 


The boundary conditions( 5,25) yield 


B 


3 



JL 1 

8a sinh(~a72) * 


(5.28a) 


Other boundary conditions (w^ = Dw^ = 0 at z = + 1/2) 
determine the constants and B 2 ,vhich are 

1 . ^3 

B.| = “ 2 ^®2 ■taoh(a/2; -i- ^) , (5,28b) 

2 R a^% (1 -k - ■“) k A 

®2 ” Ttoh-^F a T [ ;3 

x cosh(a/2) . (5.28c) 

Inserting Eqs.(5.22) and (5.27) in Eq.( 5.20b), we obtain 
2 2 

(lE +a )cos'rtz 


= _ [ c.osh(az) B 2 Z sinh(azj + B^z cosh(az) 

— cosTtz ] . ( 5,29) 




Ro - 1 ^ |) 

( it + a ) 
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Multiplying throughout with cosnz, integrating from 
z = - 1/2 to 2 = + 1/2, using Eqs. (5.28a) - (5.28c) and 
s implying we obtain 


^0 = 


,2 2.3 . 2 

(it , + a ) / a 

[(1 - k -~f)~G7(a) +|G2(a)3 


(5.30) 


where 


G^(a) 




16it^a cosh^(a/2) 

rj ^ 

^ + a ) ( sirli[a+ a) 


(5.31a) 


and 


G2( a) 


4 

a 


+ 


(-" . a") 


X coth(a/2) - cotha/ 23 < 


(5.31b) 


The expression for [ Eg. (5.30)3 , v.-hen minimized with 
respect to the *a', yields minimum value of the wavenumber 
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thereof,, critical value of Rayleigh number (7 --q^) 
in the absence of modulation. Corrections due to the 
modulation are obtained by requiring that Has ,(5,1 8a) — 

(•• 5,1 8d) and Eqs 5 ,19c )-( 5 .19d) must have solutions 
compatible with that of Egs .( 5,17a)-( 5 .17d) . The solva- 
bility criterion for Eqs *( 5 .1 8a)-(5 .1 8c) yields 

R^a^(l-k)e^ - 
-f(z,u) e ■ 

0 

/ (5.32) 

Cs> fv> CO 

where w^, 0^, are left vectors of differential operator 
of Eqs ,( 5 .17a)-( 5.17c) and ’T’, the time period of the 
modulation. The left vectors in the absence of modulation 
are given by equations 


\ 



(D^- - a^j^w^ = C (5.33a) 

-R^a^(l-k)(;5^ + ^ i (D^-a^)j^ = 0 (5.33b) 

and 

R^a^ jjw^+(D~a;j^ = 0 (5,33c) 

with boundary conditions 

w = Dw = 9,- = Dj = 0 at z = + ^ • ^ 5,33d) 

O 0 O U ^ 

We solve Eqs .( 5.33a)-( 5 .33c) by making an ansatz on 9^ ( as 
done before) > which is 
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"" ^0 • (5.34) 

i 

Inserting this in the nq.( 5,33a) and solving the resulting 
equation with boundary conditions on we find 

Wq = A^[P cosh(az} Qz sinh(az) - — — ], 

+ a ) 

(5.35) 

where 

P = - ^ tanh(a/2) (5,36a) 

and 

2 % cosh(a/2) 

Q = _ — ™ — ^ (5,36b) 

( 71: +a ) ( sinha+a) 


The evaluation of is not needed as the term .due to 
modulation in the equation for is neglected. Now fran 
Eq,(5,32), we obtain 

R^a^d- I^<W q! 0^> - ^'^o’^'o^" ^®o ^ f(z,w)^“'^w^ >= 0, 

(5.37) 

where, bars denote tim.e average over one cycle of modulation. 
As f(z,m), 0Q, and w^ are time independent, the time 
averaging of the last term of Eq.(5.37) vanishes. Therefore, 
it leads to 

= 0 , (5.38) 

which means there is no correction to the Rayleigh number 
Up to 0(e ) . To obtain the explicit forms of and 3^ upto 
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0( » f( z, o))w^( z) is expanded in Fourier series and only 

the first term is retained. Thus, 


f ( 2 , q) W|^( z) - f(z,w) V'j^( z) ^ F( a) cosnz , 

where w^(z) = — , 

and 


(5.39a) 


o 


F(w) = 


+ 1/2 

2 / f(z,«)w eosTtz dz 

- 1/2 ° 


+1/2 Y d Cosify^dC^ - z) ? 

= 2 Re f [B. coshfe^z)^B^zsinh(az) 

-1/2 sinh(Yd) ^ 2 

+ BgZ Co3h(az ) + 

R a^ (1 -k - ■^)A^ cosicz 
, 0^ s o a 

— J COSTCZ flz 


. 2 , 2.2 
(it a ) 


= 2 Re 


+ 1/2 


sinR7d)[®1 Cosh(az) -cosh^ydC-l -z)}cositz dz 


1 

nh(az) coshtydCo - z))e' 0 sitz dz 


+ 1/2 

+80 ^ Z SI 

- 1/2 

+1/2 2 1 Y, 

+ Bo z cosh(az) coshtydC^ -z/ld-osnz dz 

- 1/2 

R„a2(, -.-fAp 




1 2 

fcosh{Y=(o “•z)}cos %z dz , 

- 1/2 

(5.39b) 


Performing integrations in the above equation, we obtain 





H") = [ B, 


+ B 0 
2 w 2 


^ + 
3 3 


( 7 t ^ + a ^)^ 




4 ' 
( 5 , 40 ) 


\/vhere, 



2 yd 

® ir*(r d ) 


cosh( az) ,Sosl^d('^ - zjfcositz dz 


_.Ii_Y_d .... j- cosh{(a +Td)/ 2 > ^ __fosht(a -Yd)/ 2 > 

S'inh(^) {(a +Yd)^ +%^ } {(a -Yd)^ +%^} 

( 5 . 41 a) 



2Yd 2- ^ 
s inl^d ) da 


L ^ cosh(a 2 ) Cosh%'d(^ - z)}Cos^z dz ] 
- 1-/2 ^ 


Yd 

2 dinh('^) 


siniita+Yd)/ 2 } , sinhC( a-Yd)/ 2 ? 
C("a-i-Yd)^+ {(a -Yd)^ + Tt^ } 


4 ( a-^Y d) Co$h^(a-fY_d )/2 ) 
{(a +Yd)^'''^~^^ 


4 ( a-Y d ) Qos 1 .iCa -Y d ) / 2 ? 

{(a-Ydr -f- 


( 5 , 42 b) 



2 yd 

s inh (Y d ) 


X 



+ 1/2 

/ 

-1 /2 


cosh( az)cosh{Yd(' 2 -z)}cos 7 tz dz ] 


•^Yd r CoshCa+Yd)/2 l .. 8( a+Yd) sinrila +Yd)/ 2 } 

43 inh(^^) {(a-iYd) +'rc 3 ?(a +Yd) + ■* 3 


8 coshn a+ Y d)/2 i 

_ 2 > 2 ^ 
{(a +Yd) + ^ ^ 

with Yd “Y d 


8 (a +‘yd)^ cosh|a,+xd)/^ V terms 
f (a +Yd)^ + -r 


V 5 . 42 c) 
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and 


A 

'U 


9y d 1 \ 2 

SinltYd ) -z)lcos %z dz 


4% 


(4%^ +Y^d^) 


(5.42d) 


The time dependence in , 0^ and can be expressed as 


w-|( z,t) = W^( z) e^^ , 

(5,43a) 

0^ ( z, t) = ( z) e* , 

(5.43b) 

V / ' i^*^ 

jl( z,t) = J^( z) e , 

(5.43c) 


Then Eqs .( 5.1 8a)-( 5.1 8d) with = 0 become 

(D^-at(D^-b^)W.( z) = noa^(1-k)0i(.z^) “ 


(D^-b,)©|(z) = -VS,(z) - f(z,“) w^, 

0|(z) 

(D^-b^) J^Cz) = (D^-3^) “S~ ’ 


(5,44a) 

(5.44b) 

(5,44c) 


W,(2)= DV(,(z)=@,(z'j = D J,(z) = 0 at z = + $, 

(5.44d) 


where 


= 


and 


a im, 

a^ + imn 
a^ iw -g 


(5.45a) 

(5.45b) 

(5.45c) 
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Consistent wj.th the boundary conditions on(^^(z), we choose 

(^^(z) w A cosuz , (b,46) 

Inserting this in Eq,(5,44c), we arrive at the equation 


(D^-b^) J^(z) = - ^ 


Aj + a^)cos itz 


(5,47) 


The general solution of the above equation vath the conditions 
DJ^(z) = 0 at z = + 1/2 is 


(■^ ^+a^) tosh(b2Z ) 


Using Eqs.(5,46) and (5.47) in the Eq,( 5.44a), we have 


(5.48) 


2 2 

0 0 2 0 0 2 V (^ a )A. 

(D^-a^)(D -b"')W^(z) = R^a (l-k)A^CosTtz - R^a | — 5<; 

(n; +b 2 ) 


[ 


cosh(o22 ) 


b^ 3inh(b2/2) 
The general solution of this equation is 


+ Cosftz] , ( 5 , 49 ) 


V't(.,(z) = A^f sinh(az): C 2 oosh(az)+ D^sinh(bz) + D 2 posh(bz) 

\ 


k (52^) 


+ 


R^a (1-k) cosTtz 2 k 

o — _ _ — 27 r / 2 ' 


("rt^+a )(n^+b^) 


o s 


(n^+bj) Cb2(b2-a")(b2-b^)j) 


COS-JtZ 

(7i^-i-a^)( Ab^) 


b2/2 )} 


rl 


(5.50) 


Because of the symmetry of the system, we anticipate 

r-. = D. = 0 . : 7 (5.51) 


DU 


Now, the boundary conditions W^(z) = 0 at z = + 1/2 gives 


2 . 2 > 


C2Cosh(a/2) + D2 cosh(b/2) = 


coth( h^2) 


(TtJ+bJ) ‘■b^(bJ-a^)(bo-b^) 


"2 

(5.52 a) 


and DliV^ ( z) — 0 at z — 1/2 gives 

a C2 3inh(a/2) + b sinh(b/2) = 


Rpa (^-k) 
(•!x^+a^)(7c^+b^) 


+ 


R^a^ic (Aa^) ^ ^ 

-TZT, "ri ''2„,5"r5 


1 


s (71 +b2) 


(b2-a^)(b2-b'') {% ^+a^)( 7t^+b^) ^ 


(5«5-2b) 


Eqs,( 5.52a) and (5,52b) can be vffitten in matrix form as 
/cosh(a-2) eosh(b/2) 


( 

\ a9inh(a/2) bsinh(b/2) 


where 


cx = 


I 


+ a ) 

2 2 
(it"' + b2)b2 



I 

I 


(5.53) 


TT[^~Z^T7Tr2~"””'ir2^ 

(b2-a )(b2 - b ) 


(5,54a) 


and o 

B V 


o 9 9 

R a^ k (Tx +a ) 
o 

+ .2 . 2 " =< 


X [ —-9 — n- ' 2 ' — 2r 

(b|^(b2 - b ) 


s (n + bg) 
(it^-;-a^)(7c^ + b^) 


■ 3 . (5.54b) 


From Eq, ( 5,53) , 
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/ K / 

/ ^2 \ / Cosh(a/2) 


\ D2 


\^as inh( a/ 2) 


Cosh( b/2) 


bsinh( b/2) / 


P / 

^ / (5.55) 


which yield 


and 


b 3inh(b/2) - cosh(b/2) 

[ b3inh(b/2)cosh(a/2) - asinh( a/2)cosh( b/2).l 

( 5.56a) 


- a sinh(a/2) + cosh(a/2) 

[ bsinh( b/2)cosh( a/2)-asinh( a/2)(i osh(b/ 2)1 


(5.56b) 


Rede’findng = 


(5., 57a) 


D2 = ° 


(5.57b) 


Eq. (5.50) can be written as 


f^( z) = A^[C cosh(az) + D Cosh(bz) - cosisz + Yg Cosh(b 2 z)l', 


A.| z) 


(5.58) 


where z)= [C coeh(az) + D cos( bz)- Y , 0°= ^ 

T,= [ 1 - ic- I - ^ ^ ^ 3 (5.59b) 

^ (-n; '+a^)(^ t-b ) ^,^2 ^ 
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and 


Y 


k(it^ + a^) 

- rr — X 

^ \ 


% 


^2^ C b2-a''^)( b^-b^) s inh( b2/ 2 ) 


(5.59c) 


Inserting Eqs.C 5. 39a) , (5.45) and (5.58) in Eq.( 5.44b) 


2 2 

(^ + b^)A^ Cosicz = W^(z) + f(z,o)w^ 


(5.60) 


Multiplying the above equation throughout with cosijz and 

1 1 

integrating from z = - to z = + ^, v/e find 


= 


{tz + b -| ) 


+ 1/2 + 1/2 

/ W^(z)c'csn:z dz + J 
- 1/2 - 1/2 


f( z,:Yi>)w^Cios«z dz 


1/2 


A. 


Redefining, 
and u^(“) 


A-i </■ ( z) Oositz dz + F( u) 

^ -1/2 ^ ^ 


^0 " 


= lin 


^0 ® 




(5.61 ) 
(5.62a) 

(5.62b) 


1/2 

with K(a)) = 2 / W^(z) cositz dz , 

- 1/2 


( 5.62c) 


-1.- - .V > -'I 

Eq, (5.61) yields 


A 

X 


1 


R^a^y(ai)/(it^ + bft 
[1 - R a 


(5.53) 


To obtain next order correction to the Rayleigh number, we 


turn to Eqs ,( 5 ,1 9a)-( 5,1 9c ; , The solvability condition at 
this order of ’s ' yields 


1 

t 


T 

fN> 

/ ( ''V 

• ^ 


CO 

© 




-Fr 

-f(z,tt)e w. 


dT = 0 . 


(5.64) 


ia);5 


As — W^(z)e , the above equation after integrating 
over one cycle of modulation gives 


R^a^C 1-k 0 <WqI 0^> - R2a\<w^i3Q> - 5 <9^ I f(z,w)w^(z)> = 0 


or. 


R2a 


^ ^ ^ <0^ 1 f( Z,u) W^( z)> 

= ^ Re [ ^ ' 


W^I0^>- lc<wj j„>} 


] 


o -"O 


(5.65) 


Now using Eqs,( 5,20b) and (5.22), 


<©Qlf(z,a)) V/^(z)> = <©Jf(z,“) W^(z) 


> 


' <W (z)lf(z,o>) W,(z)> 

(Tc^+a ) 

<W^( z)I f ( z,Cj)) Wl( z) > 

(tc +3 ) 

R„a^>^.*(o)) , ^ 

2 . < ®o * '*^1^ 

(it + a ) 

R aV*(w)Ai - ^ 

— ^<© (z)lW,(z) ■>. 

( it^ + a^) ° V 


(5.66) 
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Again from Eq,(5,35), 

w = 
0 


ra 



(5.67a) 


Vi/here 


<N> 



P cosh(az) + Qz sinhXaz) - — 

(Tc^+a^)^ 


(5.67b) 


Now, Eq.(5,65) becomes 


R2a' 


= 2 Re [ 


<e^( z) I f(z,m) z) > A 


<Wq( z) J 9q( z)> 


x{1-k(l+ — )^ J.(5.68) 


<W^(z) ie^(z)> 


Inserting Eqs.(5,66) and Eq,(5,63) in the above equation, 
we obtain 


^2^ 


2 




( V^)^~1 3(^)1 ^<Qc^z)! W^(z) > 
(TT:^’-i-a^)(^^+b^) < Wq( z) IqJ, z) > 


x[l “k^1 


<WQ(z)lj^(z) > 

<V(z)ieQ(z)> )' 


(5.69) 


In most of the liquids and over a significant parameter 
range in ^He— *^He tho ratio “ is very large, oo, we will 
find correction in the limit ■g’ > >1 , In this case, 

R sinh(a/2) ( sin h a - a) 

^1 ( sinha+a}”^ ^s-^ ^-n:^ + a^)^ 32 sinh( a/2) 


(5.70a) 
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R 2 

r .Cosh_i_^£2i _ 2 Cosh(a/2)._ -i 

(sinh a+a) 4a sinh(a/2) {%^ + 

(5.70b) 


and 




Vs •* large 


{ g.lnh(a/2) ^ -id 

ITinha+a)' 32a^§inh(a/2) ^ 


+ _.l { ■ c Qsh La/2 ) _ 2vOsh(a/2)__ 

\ s i nhsV a ) ^ / ^2 , 2\5 2 


4a'^sinh( a/2) + a )' 


IT 


’5 T . rn ~<fdl • 


8a s i'nh(a/2) 


(5.70c) 


Writing 


£' 


[(1 -ic- |) G (a) + f GgCa)], 

{%^+a^r ’ 


(5.71a) 




1 


k/s •*■ large 


Now, 


“2^; 272 * s - 

(ii; + a ) 




^ [G2{.a)-G^(a)3 (5.71b) 


< ©J z) 1 W. ( z)> = ■j K - 2 


(5,72) 


where the explicit expression for ^ is given by 
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^ s • 


1 


9 (*^ + b^)b Coth( bo/2)sinh(b/2) 


. 2 ,, 2 


+b ) X ^ lcosh(b/2)3 

(fc»2-a )(b2-b^) (TC^-a^XTC^+bl 

2 2 

_ . ui /r.\ 2 r ooth( b^/2)s>inh(a/2) rj ^ 

x -aoshl a/2) ~ [ .__2 Li..,.: „ (%^+a^) 

(bj-a^)(b^-b2) 


, 1 

x{~--2 2 

(b2“a )(b2-b 


/ 2 , 2 .] 2 ., 2, 

(.-n: +a }(.-n; +-b , 


} X Cash( a/2)] cosh( b/2) 


4ir^( Aa^)(7T^+b'"^) coth(b2/2) 
b2U +b2)(b2-a^)(b2-b^) 


(5.73) 


So, from Eq,(5,30) 

I 

-»• large 

From Eq,( 5,63) , 

^0 I ,J- large 
's 


+ a^)Va^ 

[G2(a) •- G^( a) ] 


large 


(5.74) 


(5.75) 


We also have 

< W ( z)ie„( 2 ) > = - G^(a) * ( 5 . 76 ) 

0 ° 2(Tt +a ) 

Inserting Eqs.(5. 70c), (5.72), (5.75) and ( 5.76) in Eq.(5.69), 
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we arrive a t 


k 

s' 




( V^)^I Jo)) 

(?+a^r<W^ (z)lsjQ(z)> 


(5.77) 


where we have used the fact that <W (z)|sj (z)> is a term 
k o o 

of order of - and hence 


[1 -ktl - 


“gTU) - 


G^(a) 


2 H Tt^+a^) ^<W^( z ) I j z )> 


(5.78) 


The sign of correction depends upon that of the term 
<W^( z)J sj( z)> , In different range of parameter values 
this can be either positive or negative. Consequently, 
the modulation can either stabilize or destabilize the 
conduction state, 

Vi/e have seen in Chapter II that the realistic 
boundaries makes convection with zero wavenumber possible. 
It is interesting to know whether modulation make the wave- 
number finite. The correction, if there is any, will be 
because of the second order correction to the Rayleigh 
number. It is given by 


= -[( 


3 Ro 
9a„ 


3^ R. 




9 a 


s 


(5.79) 


large 

- 0 


In the limit ^ » 1 and a^ 0, this correction is proportional 


to a^. For 0, this correction vanishes. This means 

the modulation doesn't affect the convection roll, although 
it changes the onset of convection. 


S *4 Oscil latory Convection and Modulation 

At the onset of oscillatory convection velocity 
field and fluctuations in temperature field and mass current 
are time dependent (p 0), They can be expressed as 


os 

'^o 

= Re w°®( z) 

i 

6 t 

(5.80a) 

nOS 

^0 


1 Q X 


0 

e , 

(5.80b) 

OS 

'^0 

= Re J°^z) 

iS X 

0 

e , 

(5.80c) 


where is the angular frequency of oscillations. Now 

Eqs ,( 5,17a)-( 5.17c) become 


D2-a2)(D2-b2)wi°®) = , 


O ' o 

(5,81a) 


(D^-b^)6o°^' 




= . 


= (D^-a^) 


(os) 

o 


(5,81b) 

(5.81c) 


\/vhere b, b^ , are as defined in Eqs,( 2.54a)-( 2 .54c) 

(Chapter II). Eqs, (5. 81 a) -(5,81c) are same as Eqs.(2.5l)- 
( 2.53) and, hence, critical Ray leigh number the zeroth order 

t Angular frequency (S,) in non-dimensional fom 
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IS ssmG 3s "thst f 03? unmodulQ'fcoci flow, Thot isy 

2 


’ = Ros = Ss 


+ b^Va^ 


-KT- 


(1 -k~ -) G^(a,u) + f G2(a,5)} 


3 


(5.82) 


-I 

withG^(a,S) and G^CajGJ) given ^n Eqs ,( 2, 66a)-( 2,66b) [ Chapter Hi 
The terms of are given in Eqs,( 5.1 8a)X 5.1 8c) 

and can be rewritten as 


wi 




( os) 


- 

0 

^a^(1-k) 

CD 

+ 

r(0S 

) 2, .(os) 

a jX 

r(os) 

a^(l-k) 

q(os) _ 
^0 

o 

CO 






(5.83a) 

|( os) _ 
1 

1 



O 

o 

0) 

f( z,t>>) 

imt 

e 

, (5,83b) 

(D^-a 

2 - I p)4“> = 

0 . 

(5.83c) 


The solvability criterion of Eqs .( 5 ,83a)-( 5 .83c) demand 
tha t 

/ 

os) _ g( os) ^2 ^ .{ os) 
f(z.o) ei^ 


1 J (•;Us)^os)..j(^os) )/ R(»s)a2(i-u)eX' -Rr°' 


0 

(5.84) 
os) (os) 

where ’T* is the time period of modulation and » ®o 


a nd j 


^ are left vectors of the differential operator 


L = 


(D?-a2)(D2-a2-p) _r(os) 


(D^-a^-op) 

- |(D^-a2) 


R^os) g2 


0 


(5.85) 


As the convection is oscillatory at the onset, 
^^(os) i u ° 

and can also be expressed 


K os) 
wi = 

0 

Re ^ 

0 

( z) 

i w r 
e ° 

(5.86a) 

B(os) ^ 


( z) 

i a r 
e ^ , 

(5.86b) 

^os) ^ 

0 

Re ?-) 

0 

( z) 

i 

e » 

( 5 , 86c ) 


where v/ ^ and are solutions of the equations 

o ° ° 

(D^-a2)(DW) +©^0°'^ = ° 

_r(os) yos , (D2-bf)e(-) - (D^-ah 


and 


R 




0 

(5.87b) 

(5.87c) 


with boundary conditions 

= D = 00°^^ = 0 at z = + ^ (5.87d) 

It is evident from Eq.(5.84) that * do not need explicit 


expression for to determine We determine 

expression for by making an ansatz onU°^'> which 



Ao3 costcz . 


(5.88) 


Now Eq.( 5.87a) with boundary conditions on [ Eq.(5,87d) ] 
yield 



CO / \ 

M ' ' cos Ttz 

P cosh(az) + Q coshCbz) - -—5 — ^ — 5- — , 

+a )('>^^+b ) 


(5.89) 

where 


A 


P = 


os 


[b sinh(b/2) cosh(a/2) -a sinh( a/2)& osh( b/2)] 


and 


cash(b/2) 
(ix^+a^)(7c^ + b^) 


Q 


-P 


Cosh( ?'/?-) 

bc^htb/^ * 


(5.90a) 


(5.90b) 


Now 0(®) correction to the Rayleigh number due to 
modulation is evaluated by setting the system for parametric 
resonance. For this the frequency of external modulation is 
set equal to twice of that of the system. That is, 

w = 

Therefore, from Eq,(5,84), using Eqs .( 5,81 a )-( 5,81c) and 
Eqs.(5.86a)-(5.86b), we find 


(5.91 ) 



97 


y / a^(l-k)< cos^5t.c1I 


1 

T 


/ R^°=> aW“)l j(os)> to 32 „- 


1 

f 


i 

t<%y> I f( 2 ,u)> Cos^e 


T Cos 2^ T .dT= 0 
o o 


or 


a^(l-k) r(os) 32 jUs)> 

" ■5<0o°^^’ * (5.92) 


Here and are same as J^ and W respectively 

r\ r\o i ^ 


r( os) 

0 


in Eqs,(2.55) and (2.59) (Chapter II). To evaluate the 
integral at r.h.s, of Eq,(5,92), v;e expand f( 2 , 00 ) in 

Fourier cosinie series and retain only one term 


whe re 


jjUs) cosiz , 


ios) 


(5,93) 


_ 2 J* f{Zf(^)QOs%z dz (5.94a) 

-1/2 ° 

+1 /2 

= 2 •/■ CC. cosh(az) + Cg cosh(bz)+ cosh(b 2 z) + 

- 1/2 

R^°®^a^aA cosT^z Yd coshfrd(4 - z)} 

J. O °-g— — 7 ^ 1 ■ ■ . COST! z dz , 

+ .2, i,,. 1 


(Aa )(Tt +b )(•» +bp alnh^-d ) 


(5.94b) 


That is. 
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p( Os)^(j) 


[ _^;Qsh{(a+Yd[)/2} . cosh{( a-Yd)/2i 

3i«hTr-^ , i (-3.,a)2 ■; . 2} 


Txy d (ioshC(b+Yd)/2} cosh{(b-Yd)/2 ^ 

+ -« -- . — =— [ — ^ — . ■ - -2+ 

sinh(Yd/2) {(b+Yd)^+ } C (b -d)^ + ^ 


TCY d C- coshC(b^+'l(’d)/2} cosh{(b„-Yd)/2^ 

i_[ ^ +_ 4--Tr-> 


sinh( d/2) {(b^+Yd)^^-^ } C(b^-Yd)> 7t ^ 1 


We denote 


a^ a 


(it ^+a^)(ix^+b^)(TC^+b2) ’ (4T^+r^ d^) 


(5.94c) 


:( os) 


(^)| 

l»=^ 


Now we evaluate other integrals of Eq,(5,92) 




/ [? cosh(az) + Q cosh(bz) 

- 1/2 


A c os itz 
os 


+ a"^)( -nr + b ) 


1 A cositz dz 

•• no 


2 P cosh( a/ 2) 2iiQ cosh(b/2) 

" ^os ^os ^ + a^) (it^ + b^) 


(it^ + b^) 


2 (ix +a^)( '^ 4* b ) 
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.nserting Eqs ,( o ,90a)“'( 5,90b) in the previous equation, we 


arrive at 




> _ 


(5.95a) 


Further, from Eqs ,(5,89 )and(2,55) [Chapter II] by making use 
of Eqs, (5,90a) and (5.90b) ve ootain 

<wi'^)ls [H+ . . ^], (5.,5b) 


(ir^ + bS) 


v\tiere G. is given in Eq,( 2.66a) and 


71 ^ 1 

— . < rr , X 

b^( •Jt +b )(n: +b^) Cb sinh(b/2)cosh(a/2)-asinh(a/2)cosh(b/2)5 


ainh 


(bg+b) 


(b^ + b) 


sinh 


( bo+a ) 


( b^ + a ) 


sinh 


(b.-b) 


(b^ - b) 


cosh(a/2) 
sinh( bg/S) 


( b^~a ) 

•Sinh Cosh(b/2) 

—J—} ]. 

(b2 - a) sinh(b2/2) 


coth(b2/2) 

+ 2 — ^ ^ >, 

(it + bo) 


(5.96) 


Inserting Eqs, (5.93), (5,95a) and (5,95b) in 
Eq.(5.92) we find explicit expression for the first order 
correction in Rayleigh number due to modulation as 


)( os) 


k + a 




G, - 2 ^ H 


(5.97) 
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This oxpisssion is dorivod fo3? 9 03 

o* 

It IS worthwhile to look for the asymptotic limit 
( w ■* 0). In this limit, from Sq.( 5.94b) and Eqs.(2.60) - 


(2.61c) 


® b-*o 


r(os) ^2 


■ 5 5-5- (1-k~k/s) , (5.9Ba) 


Further, 

H (w) 


m -» 0 


coth(a/2 ) 


and 


(5.98b) 


G.(w)| ~ G(0) = 1 - ^ (5.98c) 

■* 0 (n +a ) (afeiriha) 


Inserting Eqs,( 5,98a) - (5,98c) in Eq,(5.97), we obtain 


R 


(os) 


(1 - k- I ) 


( j >*0 


1 -k-|}G( 0 ) - I -f4- 


(5.99) 


In case of free boundaries G(0) becomes equal to unity 

and other term with ’k/s’ does not exist. Thus 

R(os)j ~ „ i R^°®^ , which shows a destabilization. 

The result; is the same as that obtained by Agarvral et.al. for 
binary liquids with free boundaries using a Lorenz-like 
truncation of the system. However, in case of rigid 
boundaries, the explicit expression is not reported earlier. 

The correction in the limit k/ s » which is true for most 
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of the liquids and over a significant range ia is ^ 

found to be 


p,( os) 


0 )-» 0 


k 

s' 


la rge 



R 


(os) 




G(0) + 


in? 


(11 ^- 4 ^ 2 ) 


coth( a/ 2) 


which shows a destabilization 
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CHAPTER VI* 

instability in modulated couette-taylor flow 
6 * 1 Introduction 

The stability of modulated Couette— Taylor flov/ is 
a problem of long standing. The reported theoretical*^ ^ 
ano ©xparimiantal rasults do not agrse ovar the entire 
range of frequency of modulation. The full hydrodynamic 
equations for the flow have been studied analytically by 

1 T 

Hall and numerically by Riley and Laurence , While 
1 . 

Hall finds destabilization of flow due to modulation for 

all frequencies, Riley and Laurence notice, for high 

amplitudes of modulation, destabilization in lov/ frequency 

regime and stabilization in high frequency regime. For 

low amplitudes they do not notice the stabilization at 

high frequencies. The experiments of Donnelly^ show 

stabilization at all frequencies, #iile the more accurate 

experiment of Thompson shows a destabilization at low 

frequencies and a stabilization at the higher values. The 

8 

success of Chandrasekhar's technique in the analysis of 
convective instabilities in Rayleigh— Be nard flows (both 
modulated and unmodulated) and in unmodulated Couette— Taylor 
flow motivates us to investigate the stability of modulated 
Couette-Taylor flow with proper rigid boundary conditions 
using this technique . 


* Contents ofWr chapter have 


accepteo in Phys , ReV, A, 


6 .2 Hydrodynamics of the Couette-Taylor Flow 


A fluid layer is confined between two coaxial yertical 
cylinders of infinite length, the outer of which is 
fixed and the inner one is rotating vdth uniform angular 
velocity about the common axis, Havier-Stokes equations 

for yiscous incompressible fluid, in cylindrical polar 
coordinates, following Chandrasekhar's nomenclature read 
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Here Ur» and U^ are respec 


tively radial, azimuthal 


+ c nf the velocity field (U); 5^ is 
and vertical components of the v 

• d^ncitv In addition to these equations, 
the kinematic viscosity. 

have the equation of continuity. 
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( 6 . 2 ) 


The steady state for the system of Eqs ,( 6.1a)-{ 6.1c) 
described by 


is 


U 


= 0 and = V(r) 


( 6 , 3 ) 


if 


1 dp 
W Hr 


v2 


(6.4a) 


and 


V ( ^ V - (fj + i)v = 


0 . 


(5.4b) 


The expression for V can be found out by solving Eg, (6,4b) 
with proper boundary conditions on it. 

Now we consider disturbance in the steady state arxi 
assume various perturbations axismmetric and independent of 9 , 
If u^, UqjSnd u^ be the variation in steady state values of 
r, cp and z component of the velocity respectively and 6P , 
the fluctuation in pressure field, the hydrodynamic equations 
in linearized form are found to be 


3u, 

rr 

_ 2 V u _ _ 1 + 

^ r 9 P 3 r 

Vivf u^ - 
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where 

+ i 3 3^ 

1 ““ r 3r „ 2 • (6,5d) 

The ecfueiLion for conf.inui'tY fekes the "Fortii 
3 Uj. a u 

ST- + T “at = 0 • (6.6) 

The disturbance are analyzed into noriDal modes and the 
solutions of Eqs ,(6,5a)-( 6,5c) are vnritten in the form 



= u(r,t)e^^^ , 

(5.7a) 


= v(r,t) e^"^^ , 

(6.7b) 


= w(r,t) e^ , 

(6.7c) 

6P 

= p(r,t) e , 

(6.7d) 


where, k is the wave number of the disturbance in the 
axial direction,- Inserting Eqs .( 6,7a )-( 6,7d) in Eqs, (6,5a)- 
(6.6) and then eliminating w and p from resulting equations 
we arrive at following equationss 



+ ^) 
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( 6,8a) 


and 
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^ dr 
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We proceed further by considering the narrov; gap 

approximation — the separation (d) between the cylinders 

is small compared to the radius r^ of inner cylinder. In 

this approximation, -i- —) can be replaced by ^ , 

V. 

Also Q ( = ”) , in the absence of external modulation can 


be written as 

Q = 


[1 


(r - rp 

J, 


(6.9) 


Thus, the relevant hydrodynamic equations for axisyimetric 
flow in non-dimensional form using Eqs ,(6 ,8a)-( 6,8b) 
after proper rescaling are 


( ° 

- 

or 


2 

a - 


o )(L-_ _ u = (-1 V ('6.10a) 

ar 


Or 


and 


^ \ T 2 

- a - = - T a u , 


(6.10b) 


where u and v are the dimensionless radial and azimuthal 
components of the velocity fluctuation, a is the non- 
dimensional wave number of the roll. In narrow gap approxi- 
mation, there is only one variable 1= (r-r-|)/d, v\hich 

ranges from 0 to 1 . T is the dimensionless Taylor number 
defined as 


2 Q ^ d^ r 
o 




( 6 . 11 ) 


where ” = ( + r^)/2 

is the mean radius of the sys tem, v , the kinematic viscosity 
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and Xf dimensionless time paraiiKter, 

Now, vie start modulating the angular velocity of the 
inner cylinder as 

Q = Qq (1 +ecoswT ), (6*12) 


where e and o> are amplitude and frequency of the modulation 
respectively. In the steady state, the flow has only azimuthal 
component of the velocity given as 


V 


1- ^-i- eRe 


3inh{ad( 1 - 
sinh (ad) 


(6.13) 


where 



2^ 

V 


The linearized equations for the modulated flow in non- 
dimensional form are 


3 ^ 2 




a" ) u = [1 


sinh Ccd( 1-g)J - 

g+ E Re ■ r aV - ^ 
Sinh (ad) 

(6.14a) 


^ - I?’ Cccd(1-S)}ei‘^]u 

3 S 

(6.14b) 

The difference from the hydrodynamic equations for the 
modulated Rayleigh-Benard problem lies in the extra terms 
on the right hand side (rhs) of Eq,( 6.14a) as compared to rhs of 
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Eq, (4,12), The qualitative difference in the structure of 
Eq. (6,14a) from that of Eq. (4,12) persists even in the 
absence of modulation, while the modulation affects only 
one of the equations in Rayleigh-Benard system, the effect 
shows up in both equations in Couette-Taylor configuration 
and that causes a qualitative difference in the final result, 

6,3 Mathematical Analysis 

Equations (6,14a) and (6,14b) are solved perturbati- 
vely, as in Chapter IV by employing the modulation amplitude 
s as the small parameter. Expanding variables u and v 


and the control parameter T as 

2 

u = Uq +eu^ +e u^ + ...... (6,15a) 

V = +ev^ +e^V2 + (6,15b) 

T = Tq +eT^ + (6.15c) 


and then inserting in Eqs,(6,14a) and (6,14b) and equating 

/ 2v 

the identical powers of e , we have up to 0(e ) 


(D^ - a^ - 

|_)(d 2 - a^) = (1 -S)v^ , 

(6,16a) 

(D^ ~ - 

% = -vS: 

(6,16b) 

(D^ - - 

g~)(D^-a^)ii, = (l-^)v^ + Re g(§ )v^ 

i(j>T 

e , 


(6.17a) 
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,2 2 3 


“l - T, a" - T„a"u^Re fCSje^”*, 

(6.17b) 


±m 


,2 2 9 9 9 


»)(D"-a")u2 = (1-^)V2 + Re g(g) v^ , 


(6.18a) 


,2 _2 3 


(D- - a- - = -T^a ~ T^a^u^ 


where 


g(r/ = 


2 ■ .-v iWT 

- T^a Re fi£) e *, (6.18b) 


s inhCQd( 1 -§.) } 

Vinh (id) 


" K ^ i caTdT cosh{ad(1-4)} 


(6,19a) 

(6.19b) 


and D, 


d 

T1 


(6.19c) 


The boundary conditions require that all three ccanponents 
of the velocity must vanish at the rigid walls. This 
implies that 

^ = v ^ ~ 6 at§=0 and 1 , (6,20) 

We begin' with the' '0(}1 ) equations, Eqs.(6,16a) and 

(6,16b) which represent unmodulated flow. The boundaiy 
conditions at this order one 

du . ■ ■ ■ . 

V ~ Uq — ~ 0 at ^ = 0 and 1 , (6.21 ) 

Here, the only requirement on the azimuthal component (’^q) 
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of the velocity is vanishing at the vralls. So, we proceed, 

as before, by expanding v in Fourier series and truncating 

o 

the expansion at first term. This leads to making an 
ansatz 

^ sin 7!;^ , (6.22) 

where the undetermined amplitude is chosen to be unity. 
Using Eq,(6,22) in Eq,(6,16a) and noting that there is no 
time dependence in unmodulated flow, we find 

(D^ - a^)^ u^ = (1-^) sirnil. (6.23) 

The general solution of the above equation is 


Uq(I ) = Aq sinh(al) + B* cosh(aS) + CQ'fsinh( a ^ 

(l-S)sin('rt^) 


4” D S cosh( B. 4- ^ 


( ) 


(6.24) 


where A^, B ', and are determined by the boundary 
conditions on u^. They are found to be 


A* = 
^o 


47t( 1 + c osh a) 




Tt a 


( sinha-a)( + a ) 


72 “; 272 , . .2 ^ * 

(tt: + a } (sinh a - a ) 

(6.25a) 


B, 


4 Ti; 

’TT~^ 

{% + a ) 


(6.25b) 
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% 



~^rT 

+ a ) 


(sinh a cosh a - a) 
^ 1 . 

{ s inh a - a ) 


(6.25c) 


and 



4n;a( 1 + cosh a) 
(sinh — a)('J^ +a 



Tisinh^a 

“7 SrSTTT^ 5. * 

+ a ) (sinh a-a ) 


(6.25d) 


Inserting this solution Eq,(6,24) in Eq,(6.16b), multiplying 
by sin('TC^), integrating fron § = 0 to 1=1 and then using 
Eqs,(6,25a) -- (6,25d), we find that 


T 


o 


2(71 + 


2.3 
a ) 


.-1,6a.'^_^. e osJi^ (,a/.2 ) j 

{Tz+a^j{_a + sinha) 


(6.26) 


The critical Taylor number for the onset of instability in 
the absence of modulation is determined by minimizing T^ 
with respect to a. The result, as stated in Chapter X, 
is within ^% of the exact numerical ansv^r. The one— mode 
truncation is thus found to be a very good approximation 
and can be used for studying the modulated system, 

Nov 7 , before proceeding further vre determine the left 
vectors satisfying the equations 


(D^-a^)^ u^(g) 




(6.27a) 


(D^-a^)vQ( ^ ) = (1-f)u^(§) 


(6,27b) 


du. 


with the boundary condit-lons Vq=Uq - - 0 at ^ 0 and k 1 
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Again making the ansatz = sinitg, Uq(§) is found frcan 
Eg ,( 6,27a), wiiich is 


( D^-a^) ^ u^( ^ ) = -1’ sin ti g . 


(6.28) 


The solution is given by 


Uo(£) 


= A^sinh( a I) + cosh(aS) + sinh(aS) + 

+ s cosh(a s) - -*-75 — ry^ sin iis , 

° (ic^+an^ 

(6.29) 


The boundary conditionson u^( ^) yield 


ro 

B = 
o 


D. 


c 


= 0 


To" 

272' 

( 71 +a ) 


sinh a 


(a + sinha) 


(6,30a) 

(6.30b) 


and 


a Aq = 


To"^ 

(71 + a ) 


~ D, 


(6.30c) 


The solvability condition for Eqs,(6,17a) and (6,17b) 
yields T^ . The condition is expressed as 


(^0 ^0) 


Re g(^ )V q e 


rmr 


\ = 




0 , (6,31) 


where the bar denotes time averaging over the time period (T) 
of modulation. The above equation is rewritten as 
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F = 


“d 


^ ^ ^ r r i^'ha'd ) ^ os hCa d(l - £) }] X [A^ sinh(a§) + B^cosh(ag) 
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31 Q Q O O 
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;(;:7d7'«T^' 


X 
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^0 I- sinh(ad) 
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T ^ .V, f g( ^)v e 

/ (u„ v^) ^ ° 

O 


o 




dT = 0 


or, 


T 


UqI g( S)Vq> “ ^ •'‘q coscot 


T 

/ cosorc d'C = 0 

o 


which leads to T^ =0, 

To solve the 0(e) equations, vje first introduce the 
Fourier expansion of the functions g(^)v^ end f(^)u^ as 


and 


where 


and 


gCDv^ = G sinCitg) (6.32a) 

f( l)u^ = F sin( , (6.32b) 
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r cosh(o^d) - 1 I 
^ a d sinHfcd ) 


(6.33a) 


F 


1 

2 / 
o 


f( g) sin(ic^ ) dg . 


(6.33b) 


Inserting expressions for f(^) and Uq(^) from Eqs,(6,19b) 
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For 0) 0 


F( 0 [)= 0 ) 
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(6.35) 


while for 0 ) -► 00 


F(o>) 


(ad)' 


(6.36) 


The solutions of Eqs,( 6.17a) and (6.17b) proceed by choosing 
and as follows s 


'1 = 


(6.37a) 


= u^(s) e' 


(6.37b) 


with the ansatz v^( ^) = sin(7c§ ) 


(6.38) 


Inserting the above expressions in Eq.( 6.17a) and using 
Eq,( 6.22a), we have 

(D^-b^)(D^-a^) u^(^ ) = A* (l-§) sin(7ig) + G(w ) sin^t^ , 


(6.39) 


AAihere 


2 2 
b = a + iw 


(6.40) 
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The general solution of Eg, (6.39) is 


= P cosh(al) + Q sinh(al) + R cosh(bl) + s sinh(b ^ 


sin(7rg ) 
+ b^) 


2ftA*(2'n;^+a^+b^)cos(7t 
(tc + a ) (h +b ) 


G . sin( Tt 

(7t^+a^)(7t^ + b^) 


(6,41) 


The boundary conditions on u^(§) are 


u^(g) ^ 


The use of 


du^( t) 
"d t ' ' 


= 0 at § = 0 and 1 


u^(l) = 0 at ^=0 and §= 1 give 


P + R = 


2Af(2-n:^ + a^+b^) 


o "J 9 9 o 9 

(Tc^+a ) (T^^+b )'^ 


(6.42a) 


and Pcosha + QSinha + Rcoshb -i- S sinhb 

"^9 2. 9 

~ 27 i;A^ {0% -r a + b ) 


while that of 

^) 


{%^ + a^)(ic^ + b^) 


0 at ^ = 0 and t = 1 yield 


(6.42b) 


aQ + bS = 


~Tu(A^ + G) 

— IS ' - ' -g — 2T 

(it^+a^)(Tf+b ) 


(6.42c) 


aFslrha+ a Q cosh a + b R sinhb + bS coshb 


( 1 C +a ) ( ic +b ) 
(6.42d) 
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we define 


X 


* 2 2 9 

2itA^ ( 2 % + a + b ) 

' ■ 7^2 — : 7:2 » 

(tc +a ) (II + b } 


(6,43a) 


Y = 


n (a" + G) 


(Ti^+a^)( 7? + 


(6.43b) 


Z = 


It G 

TT — n- rj — rj— 

(iX +a)(7I +b) 


(6,43c) 


2 2 

and A = ( a +b )sinha sinhb + 2ab - 2ab cosha cosh b. 


(6.43d) 


Now solving Eqs .( 6,42a)~( 6,42d) , we find 


pA - X [ b^ sinha sinhb + ab ( 1 + coshb)(l - cosha)] 

+ Y (a cosh a sinhb - b sinha coshb) 

+ Z (a sinh b - b sinha ) , (6.44a) 


QA = X [ab sinha (iHcoshb) - b^ sinhb (1 + cosha)] 

+ Y(b cosha coshJo - a sinha sinhb - b) 

+ Z( b cosh a — bcoshb) , (6,44b) 

R A = X[a^ sinha sinhb + ab (1+cosha)(1 - coshb)] 

+ Y(b sinha coshb - a cosha sinhb) 

+ Z(b sinha - a sinhb), (6.44c ) 
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SA = X[ ab sinhb+ ab sinh b cosha - a sinha Coshb 

— a^sinha ] 

+ Y( a cosha coshb - b sinha sinhb ~ a) 


+ Z( a coshb — a cosha) , 


( 6.44d) 


Now, using Eqs.( 6.32b) in Eg. (6.17b) 


(D^-b^) v-j( ^ = -T^a^ u^(S) - T^a^ F sin(^^) 


or, 


(it; + b^) sinCii; s) 




^a [ u^( I ) + F sin( TtS)] , 


Multiplying throughout with sinCiil) and integrating over t- 0 
to g = 1 
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= ^'[ / u^(^) sin(itg)d^ + F/2] 
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Inserting for P,Q,R,S and then simplifying 
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(6.45) 


So, 
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T a F 
o 


( 1 + I.X) 


(It ^ + b") 
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TTi^' 


1 - 


2(11 ^-i-b^) 




(6.46) 


where 
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(TU 4-q )(^ +b ) 


X 


{ asinhb(l + cosha) - b sinha (1 + coshb)}] (6*47) 

2 

We now turn to the 0(e) equations. The solvability condition 
at this order using Eqs..(6,18a) and (6,18b) with = 0 
yields 



(6.49) 

The explicit form for the fractional correction because of 
modulation is f ound to be 
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^2 
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= Re [ ^ 


2 2 
Tq a"" I FI 


- 21 


^ 1 

2 


f'"" G j^Cir^ + b^) 


< u (g ) lu (^> - 2(Tt'‘^+b^)l 


r2,K2K 


o 


+ 




{ 


{% - + a) 


F - G/ 

T^a”-^- 2(7c'^ + b'^) 


1 H*] 


where 


H = / f(^) “u ( ^ 

* o 

and H is its complex conjugate. 


f 


(6.50) 

(6.51) 


^ Results .in Asymptotic Limits 

It is worthwhile to examine the expression in two 
limits of low and high frequencies. 

In the low frequency regime ( a> -* 0), 


f(g ) -*• 1 and g( S ) -*• 1- ^ . 


Now, 

lim 
m “* 0 


< u^f g( § ) 


■2^ 

a <v. 


v^( S)> 

u > 
o 


lim 
0) -► 0 


/ u^(§ ) v^(^) dg 
o 

a^ / Uq(^) d§ 

o 


= lim 

(j)-* 0 

- - lim 
m-* 0 


/ (D^-a^)vQ(g)v^(S )dl 

a^ / VQ(g) ujS) dg 
o 

* 


(6.52) 
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Further 

lim 

« -* 0 


To<vj f(S ) u^(£ )> 
c o 


= lim 
m '*• 0 


T ^ V ( I) u.(g ) d I 
o o ° ' 

u fS) 


1 


- lim 
0) 0 


/ (D^-b^)v^( ^)Uq( S)d5+TQ a^/ v^C^u^l^dS 

/ v^( g) ) d^ 

o 


lim r -T + ( 7^ . ^ + fat A* T 1 

L T„ ^^2 ^2) A, • 


(6.53) 


m -* 0 


Inserting these equations Eq,(6,52) and (6,53) in Eq,(6,49) 
T 2 (t^= 0) = lim [1 ~ 


im A. 

^ - 4] . C6.54) 


0) -*• 0 


^ 2 ^ 2s 

(it + a ) 


-'■i. 

Using expression for A^' Eq,(6,46) , vie arrive at the limits 

and for a cj 3,12, we find 


0 ) 


ro -0.10 , 


(6,55) 


which is in good agreement v/ith Hall’s value of -0,07 , 
In high frequency ( m ^ ) regime, 


G/F 


1 


[ 


3a 

T 


(ic^ + a^)^ ""(sinh"^ a-a) 


t s inha cosha— a 1 — 


4 a^ (l + cosha)(sinha) 4a sinha( a cosha-sinha) 

2 . "‘I-T / 2 . 2. 


- 1 3 


(it ^ + ai 


(n^ + a ) 


is a constant 
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< 7 ^ ■* thus has the same form as the corres- 

ponding (A^/Aq) in Eq.(4.4l) for the Rayieigh-Benard problem. 
The last term in Eg, (6.49) dominates at high frequencies 

of modulation because and u^ of the first term and u and 

O i o 

of the second have similar structures and G is small than 
F by a factor This suggests that the threshold shifts 

in the Couett'^-Taylor and Rayleigh-Benard problems have 
identical structures in the large frequency limit. For 
a 21 3,12, we find 



1 .2 

5 


(6.56) 


1 

This is in disagreement v/ith Hall's analytic calculation 

which shows T^/T . The difference is not only in 

Z 0 

the asymptotic power law, but also in sign. It is the latter 
which shows the qualitative change in stability pattern of the 
modulated Couette— Taylor flow. 


6.5 Discussions 

2 

In low frequency regime {q ^ 50 in units of v /d ) 

the result shows a clear destabilization^ and quantitative 

1 

agreement with Hall is excellent. In high frequency regime 
we find stabilization -- — a qualitatively different result 
from that of Hall, Even the magnitude of the effect has 
different asymptotic dependence on the frequency of modula- 
tion, This is clear from an observation of the hydrodynamic 
equations [Eqs.( 6.14a) and (6,14b)] for modulated flow. The 
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0(e) term in Eq.( 6,14a) is smaller than the 0(e) term in 
Eq,(6.14b) by a factor ad. At high frequencies a d •* oo , 
and thus the modulation effect \vill be present in v- equation 
alone and the structure would be the same as that of corres- 
ponding Rayleigh— Benard problem. And, therefore, stabiliza- 
tion is not uneypected. The same qualitative result is 

9 

obtained by a Lorenz— like truncated system for CT flow, 

Riley and Laurence’s numerical procedure is not 
sensitive enough to find the shift in Taylor number at such 
high frequencies for low e , However, for s>1, they do 
observe a stabilization at higher frequencies — a result 
in qualitative agreement vdth our work. The sudden dis- 
appearance of this effect from their result, as e is lowered, 
is probably due to the smallness of the magnitude of the 
effect and the consequent numerical difficulties. As for 

the experimental studies, Donnelly reports stabilization 

7 

at all frequencies. The later experiments of Thompson shov; 
a low-frequency destabilization but a high frequency 
stabilization. 

Our approach is able to produce the results in the 
entire range of the frequency of modulation for small e , 

Thus it also provides an alternative approach to that of 

1 

Hall for analyzing modulated flow. 
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APPENDIX 


STROBOSCOPIC MAPS AND CHAOS IN MODULATED BINARY LIQUIDS 

A,1 Introduction 

1 -4 

A lovj-oTder truncation of hydrodynamic systems 
facilitates the handling of non-linearity, while preser- 
ving the correspondence with actual hydrodynamic equations 

3 4 4 

in case of double-diffusive systems ’ . Agarwal et,al, 

carried out a Lorenz-like truncation for a system of binary 

liquids and reproduced hydrodynamic results regarding the 

convection threshold® , Here, we construct a stroboscopic 
5 

map associated with a modulated binary liquid near the 

onset of oscillatory convection. In a suitable coordinate 

system the only relevant direction of the map is isolated. 

The existence of a return map in this direction shows the 

possibility of a the period-doubling route to chaos . 

The threshold for the onset of chaos in the experin^ntally 

3 4 

interesting system of He— He is also ccxnputed. For 
simplicity, we work with idealized boundary conditions in 
this section, 

A, 2 Hydrodynamic Equations 

Lorenz-like truncation for a modulated binary liquid 

yield the following equations; 


X 


a (-x+y+u) , 


(A.1a) 
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• 

y = 

-xz + r ( 1 

+ A cosmt) X - y - u. 

(A. 1b) 

♦ 

Z 

xy - bz - 

ti^s 

(A. 1c) 

• 

U = 

-XV - kr ( 1 

+ A coswt) X - su -ksy. 

(A.1d) 

* 

V = 

ux - bsv - 

kbsz , 

(A.ie) 


The five modes represent velocity potential, circulating 
modes of temperature and concentration, and temperature and 
concentration modes of connective transfer of heat and mass. 

4 

Here, r is 2 the reduced Rayleigh number with ^ — I 

2 \ 

^ - )^T » where Cp is the specific heat and x , the 

susceptibilityiA and ware amplitude and frequency of modu- 
lation respectively! b = 8/3, All other symbols have the 
same meaning is defined in Chapter II, 


We now carry out a transformation of variables in 
Eqs,(A,1a) - (A1 .e), which is designed to make the z and v 
equations in hlinearized form version of Eqs .( A,1a)-(A1 ,e) 
independent of each other. The transformation needs 



where 

cos^e =[1 + /2 


(A.2) 


(A.3) 
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and 




(A.4) 


If we define a vector X vdth components (X^ > X^, X^, X^, 
Xg) = (x, y, u, z, v) then the system of Eqsi(A,1a) - 
(A1 ,e) can be revinrittoo as 


X 


^ = cr (-x^ + X^ + X3), 


where 


X . 


X . 


a 


and 


X 


1,2 


(A.Sa) 


Xo = 


X , 


r(1+A cosmt)X^— X^ 


-kr (1+A cosmt) X. 




-- X3~X^(X4Cos0 - X^sinG) 

(A.5b) 


ksXg - SX3 


- a X^(X^sin9 + X^cosG), 


-X^X^ + X^X^, 
• 2^5 


x . x . 

-X X 4- ■■■ .1. ..3. 

+ rv t 


= k'/’b/p. 


~ i[b(l+s)+^3 


(A.5c) 

(A,5d) 

(A.Se) 


(A.6) 


A, 3 Perturbation Theory 

2 

By setting r = (1 )» where £<< 1 and 

the threshold for the onset of oscillatory convection 
(rnc<ro> "the threshold stationary convection) , vie carry 
out a perturbative solution for X^ , X^, X3, X^ and X^, 


The 
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system of Eqs.(A,5a) - (A.5e) can be written in compact form 
as 




A coswt(6.2 - k i = 0, (A.?) 


where 


/ 


WE 


+ a 



-a 

3 

sk 

0 

0 



9 a. 
^ + s 

0 

0 


0 

0 

0 

3 . 

0 



(A. 8) 


is the vector of nonlinear terms having the form 


N. 

1 


/ 


1 




-X^jX^cosG + X^Xg sinG 


-X^X^ainG + X-|Xg cosG 


X,X2 


a 

Cs> 


(A.9) 


and S . . is Kronecker delta symbol, we work with ~ 

the frequency for parmetric resonance and to apply perturbation 

theory consistently we set A = A £ • 
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We begin by noting that the amplitudes will have 


the form 


X = ew(t) X 


e ° + + ... , (A.10) 


where is found to be 


.(o) _ 


^os *■ / k) 


1 - ( P- s/k) - i 0)^ 


(im^/o) + 1 _ 


(A. 11) 


The perturbation theory provides W(t) and higher order 
terms. The bilinear nonlinearity allows to be expanded 


Ni = + .... 


(A. 12) 


Now considering terms up to 0(e ), we find that 


= 0 upto the same order in e and 


= 4 °' 


(A. 13) 


The solution is 
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( 1 ) 


0 
0 

/ \ iw t 0 

X^°)e ° „ 

(_2 )„2 ^ 

+2i“o 

/ \ 2i 0) t 

e ° 2 

(-2 —) ^ + 

x;2 + 2ia.„ 2 


2 Re X 


(o) 


2 Re X 
a X 


(o) 


IW!' 


I w r 


(A.14) 


We now turn to the 0( s ) term in Eq, (A,7) to find a relation 

ioj t 

for determining W(t). Equating the coefficients of e ° to 
zero, we find a condition on W in the form of the differen- 
tial equation 


H = (a^ + ia2)W - + iPg) WiW( ^ 

r * 

+ ■^ ( i 0^) ^ t 


(A. 15) 


where 


r = e t . 


(A.l6) 


a 


^ + i “2 and + ip^j are obtained from the scalar products 

(A.17) 


and 


+ ia^ =<x[°^)M > /< x[°^| X > 


P + ip2 = < ^ ^ ^ ^ ^ ^ 


(A.l8) 


where 


X 


(o) 


= ( 1 , 


( i^n ~ - . 

2-. K k y , Y , 0, 0), 


os 


(A.l9) 


133 


M = 


■os 



and 


N = 


N^CosG + N^SinG 

a (N^CosG - N^CosG) 
0 
0 


with 


S/U3^ ^)(0-+ jCO^) 
s( 1 + p,^) 4- io) 


N. 


X 


(o) 


+ 2i<0^ 


Re X 


(o) 


2 


and 


N, 


X 


(o) 


^ 2 + 21“o 


Re 


J\ o) 

yvo 


+ 2 


(A. 20) 


(A. 21) 


(A.22) 


(A. 23) 


(A. 24) 
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■'^*4* Stroboscopic Map 

The concept of discrete map has proven itself 
worthiAhile in analytic study of chaotic behaviour in 

5 

dynamical systems. Construction of a stroboscopic map 
is a powerful instrument to study non-autonoraous systems 
which do not admit multiple time perturbation scheme, 

A stroboscopic map in an n-dimensional non-atrtonomous 
dynamical system is a discrete mapping from an n-dimensional 
space to other obtained by considering the flow at times 
separated by regular intervals, V^e note from Eg, (A, 10) 
that the form of limit cycle is specified by the amplitude 
W, which is complex. Hence fore < <1 , we can reduce the 
5x5 flow system to a 2x2 map. 

We begin by taking the form of W( r) as 

W(t) = ^ (A, 25) 

where v^x), radius of the limit cycle, and^C"^) are real. 
Inserting Eq. (A, 25) in Eq.(A.1 5) and equating real and 
imaginary parts separately, we have 

w = w-P^ w^ + 6 w cos2(]S - 9) (A.26) 

and 

jS = a2 (^2 Ssin2 (0 ~9 ) (A, 27) 

who re 

26 = a [c^ + a 2 ^ , (A,28) 


and 


I 


Again the recursion relation for w is obtained by inter- 
polating Eqs,(A,3:0) and (A. 34) using (A.S?). It is found 


to be 


w , . = w„ e 

n+1 n 


(a^+6)T 


° [1 + Aw^] ’''2 cos(-^. 


COS c-^.) y 


,[ 1 + (^) ] . 


{A.38) 


There are only two independent components of X up to 0(£ ) , 
We construct tv/o variables as 


^n ^n » 




(A. 39) 


Using the recursion formulae (A,36) and (A.38) after 
rotating the coordinate axes by'!', 'm find 


a.T 

e’ ° 

1+(Pn-^-qnK] 


0 W cos§^ sinf^ 


-s.in^^ 6os^_ 

V n n 


where 


^ [In ( 1 + p^A)] 


(A^O) 


(a, 41 ) 


After large number of iterations the system (A,40) reduces 


to one dimensional return map 


Z 27 ». 1 / 2 

[1 + Pn^J 


(A.42) 
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Thus, we find that for modulated binary liquid the strobos- 
copic map reduces to a one— dime nsional return map in the 
0(e) analysis. The map is easily seen to have a quadratic 
top and does show a series of period doublinq bifurcations 
to chaos. 


A, 5 Numerical Results and Conclusions 

As an example, we consider the typical values, 
cr = 1 , s = 0.05, k = 0.04 and p ~ 0.01 for %e-%e mixture 
at a temperature of about 10 mK from the superfluid 
transition point. The linear stability analysis yields 
r = 1.12, r = 6.25, so that the onset of convection will 
be oscillatory. The frequency cOq of the oscillation is 
0,14, We set r = r^g (1 + ^ ), where s = 0.01, that is we 
are slightly above the threshold for oscillatory convection. 
The radius of the limit cycle 21 ^ constants 

p^andp2 -0.065, and -0.26 

respectively. The stroboscopic map of Eq, (A, 41) now turns 
out to be 


Pn+1 



cos?-^ In (1+0.4) p^) 


{1 + 0.4 } 


2 , 1/2 


(A,43) 


where 



1 .2 e 


0.36 


(A, 44) 


This map has a quadratic max at Pj^' - Pq 3.5 and shows a 
series of period doubling bifurcations. The limit point 
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in terms of A [ the amplitude of the modulation in 
Eqs,(A.5b) and (A.5c)3 turns out to be A ^ 0.2. 

Numerical integration of the dynamical system for 
the above parameters shows a sequence of period doubling 
bifurcations after some precursor transition where the 
limit cycle of frequency changes shapes. These 

precursor^ transitions cannot be handled in terms of the 
stroboscopic map introduced by us in Sec,(A*4), The 
threshold for chaos in the numerical experiment rs A ^ 0.38 
in qualitative agreement -with A ^ 0.2 obtained from the 

0(e ) map. 
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